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ABSTRACT 

We present Herschel/HIFI observations of the three ground state transitions of H2O (556, 1669 
and 1113 GHz) and H3, 8 (547, 1655 and 1101 GHz) - as well as the first few excited transitions 
of H2O (987, 752 and 1661 GHz) - towards six high-mass star- forming regions, obtained as part 
of the PRISMAS (PRobing InterStellar Molecules with Absorption line Studies) Guaranteed Time 
Key Program. Water vapor associated with the translucent clouds in Galactic arms is detected in 
absorption along every line of sight in all the ground state transitions. The continuum sources all 
exhibit broad water features in emission in the excited and ground state transitions. Strong absorption 
features associated with the source are also observed at all frequencies except 752 GHz. We model 
the background continuum and line emission to infer the optical depth of each translucent cloud along 
the lines of sight. We derive the column density of H2O or H^O for the lower energy level of each 
transition observed. The total column density of water in translucent clouds is usually about a few 
10 13 cm -2 . We find that the abundance of water relative to hydrogen nuclei is 1 x 10~ 8 in agreement 
with models for oxygen chemistry in which high cosmic ray ionization rates are assumed. Relative to 
molecular hydrogen, the abundance of water is remarkably constant through the Galactic plane with 
X(H2 0) = 5 x 10~ 8 , which makes water a good traced of H2 in translucent clouds. Observations of 
the excited transitions of H2O enable us to constrain the abundance of water in excited levels to be 
at most 15%, implying that the excitation temperature, T ex , in the ground state transitions is below 
10 K. Further analysis of the column densities derived from the two ortho ground state transitions 
indicates that T ex ~ 5 K and that the density n(H2) in the translucent clouds is below 10 4 cm -3 . We 
derive the water ortho-to-para ratio for each absorption feature along the line of sight and find that 
most of the clouds show ratios consistent with the value of 3 expected in thermodynamic equilibrium 
in the high temperature limit. However, two clouds with large column densities exhibit a ratio that 
is significantly below 3. This may argue that the history of water molecules includes a cold phase, 
either when the molecules were formed on cold grains in the well-shielded, low temperature regions 
of the clouds, or when they later become at least partially thermalized with the cold gas (~ 25 K) in 
those regions; evidently, they have not yet fully thermalized with the warmer (~ 50 K) translucent 
portions of the clouds. 



1. INTRODUCTION 

The ratio of the ortho to para spin modifications of a 
molecule containing identical atoms is a valuable tracer of 
its chemical history. The original formation of the species 
in question, as well as its subsequent history, affect the 
relative abundance of the molecules having spins parallel 
(ortho) and antiparallel (para). The different spin modi- 
fications are treated as almost entirely separate species as 
radiative transitions are generally highly forbidden, and 
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collisions that convert one spin modification to another 
are very slow, with their rate reflecting the density and 
temperature experienced by the molecules throughout 
their history. The observed ortho-to-para ratio (OPR) 
is thus sensitive to physical conditions throughout the 
lives of the observed molecules. 

A variety of effects complicate using the OPR to probe 
the history of molecules in dense clouds. Excitation 
and radiative transfer effects make it difficult to com- 
pare accurately the column densities of the different spin 
modifications. A complex chemical network with numer- 
ous formation and destruction pathways make it diffi- 
cult to disentangle to what other molecules' OPR tha t 
of a given species is related (e.g. iKahane et all 11984] ). 
The physics of spin exchange processes for binary col- 
lisions, but even more so for surface interactions (e.g 
iBuntkowskv et ahl 120061) . is very complex and not fully 
understood. Despite these challenges, the OPR has been 
observed using molecules in a variety of conditions and 
the observed ratios used to model current and past con- 
ditions. For example, the rotational transitions of molec- 
ular hydrogen are re adily observed in shock-heated gas. 
iNeufeld et all (jl998l ) studied H 2 in HH54, and derived 
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shock conditi ons from the ob s erved OPR. Subsequent 
modeling by IWilgcnbu s et al.l (|2000D indicated that a 
complex set of shocks was required and that the solu- 
tion was not unique. Much bett er data was obtained 
subsequently (Ne ufeld et alJl2006[ ). but interpreting the 
significance and time dependence of the OPR remains a 
challenge. Studying molecular absorption in interstellar 
clouds is a potentially powerful probe of the OPR, as 
well as of the physical and chemical conditions in these 
regions. As discussed at length elsewhere, the relatively 
low density in these regions makes the collision rate mod- 
est, and thus for many species it is a good approximation 
to assume that only the molecular ground state has ap- 
preciable population. Since the background sources are 
generally of small angular size, at velocities for which the 
absorption may be optically thick, the line should have 
a flat-bottomed appearance, which is a valuable confir- 
mation that the continuum level has been determined 
accurately, and helps define the velocities for which the 
emission is optically thin. 

Water, like H2, is a molecule with an ortho and para 
form. Water vapor has been widely observed in dif- 
fuse and translucent interstellar clouds, by means of ab- 
sorption line spectroscopy towards bright submillimc- 
ter continuum sources. Observations ca rried out us- 
ing t he Infrared Space Observatory (ISO, Kessle ret al.l 
19961). Submillimete r Wave Astronomy Satellite (SWAS 
Melnick et alJ[l 999) and the Herschel Space Observatory 
(jPilbratt et al.ll2010T ) have revealed typic al water abun- 
danc e s of a few 10~ 8 relative to H 2 ( e.g. iNeufeld et al.l 
12001 1201tt iSonnentrucker et al.ll2010h . However, signif- 
icantly higher water abundances, of a few 10 -7 , have 
been reported for diffuse molecul ar and translucent 
clouds in the Galactic Center region (INeufeld et al] l200"fT: 
iCernicharo et al.ll2006t iSonnentrucker et all 120 nt h 

In the diffuse interstellar medium, the destruction of 
water is dominated by photodissociation. Three sep- 
arate water f ormation mechan i sms have b een posited 
to date (e.g. iHollenbach et al.l 120091 120121 and refer- 
ences therein). Much details abo ut water in the in- 
terste llar medium can be found in Ivan Dishocck et al.l 
(| 2 1 lh in their presentation of the early results of 
the Herschel WISH key program, mainly dedicated to 
the study of water in star forming regions and proto- 
stars. First, water can be produced by means of a 
gas-phase ion-neutral chemistry that is driven by cos- 
mic ray ionization. The chemistry is initiated by cos- 
mic ray ionization of H or H2, followed by the forma- 
tion of OH+ via the sequence H+(H 2 , H)H+(0, H 2 )OH+ 
or H+(0,H)0+(H 2 ,H)OH+; the OH+ ions thereby pro- 
duced undergo a series of hydrogen atom abstraction 
reactions, leading to HjO+j which ultimately produces 
water vapor as a dissociative recombination product 
OH+(H2,H)H20 + (H2,H)H 3 + (e,H)H20. The fraction 
of H30 + recombinations that lead to H2O has been mea- 
sured in several laboratory experiments, the most recent 
of which - perfor med in ion-storage rings - suggest a 
value - 17-26% (|Neau et alJl2000Tl Jensen et alj 2000). 



A second formation mechanism, the importa nce of which 
was emphasized by IHollenbach et al.l (|2009f) , is the hy- 
drogenation of oxygen atoms on grain surfaces, followed 
by the release of water ice by photodesorption. Finally, 
enhancements in the production of water are predicted 



to occur in regions heated by shocks or turbulent dissi- 
pation to temperatures above ~ 200 K by a sequence of 
two neutral-neutral reactions, both possessing activation 
e nergy barriers: OfHz, H)OH(H 2 , H)H 2 0. 

INeufeld et al.1 (|200Q[ ) observed the ground-state tran- 
sitions of ortho-H^O and ortho-H^O in absorption to- 
wards Sgr B2, while INeufeld et aTT(|2002l) observed the 
H^ 6 isotopologue towards W51. iPlume et all (12004ft 
carrie d out a similar study towards W49A. iGerin et HH 
(2010) have observed CH in absorption towards W31(C), 
W49(N), and W51, obtaining information on the phys- 
ical characterist ics of the line of sight absorbing clouds. 
iLis et al.l (|2010D analyzed H 2 Q and H^ 8 Q data from the 
HIFI instrument (|de Graauw et a l. 2010) on Herschel for 
lines of sight towards Sgr B2(M) and W31(C), to de- 
termine the OPR ratio in different velocity components. 
Most of the results were consistent with the statistical 
high-temperature limiting ratio of 3. However, for ve- 
locities corresponding to the expanding molecular ring, 
this study found OPR = 2.35 ± 0.35, which could imply 
equilibration at a lower temperature. 

In t he present paper we extend the work of ILis et al.1 
(poTol ) to the remaining PRISMAS lines of sight, towards 
DR21(OH), G34.3+0.1, W28(A), W33(A), W49(N), and 
W51, which have been utilized as background sources 
previously, but only now with the capability of the HIFI 
instrument on Herschel, can the ground state transitions 
of ortho- and para-H20 be observed with the high sig- 
nal to noise ratio required to determine the OPR in the 
interstellar clouds along these lines of sight. In S}2] we 
present the observations. In <JJ] we discuss the difficult 
problem of establishing the proper background spectrum 
to use in determining the absorption optical depths. The 
reader not interested in the details of the modeling can 
skip this section. In Sj4] we discuss our results for both 
the abundance of water and the OPR in these translucent 
clouds. We summarize our results in <J5] 

2. OBSERVATIONS 

We present here all the HIFI (|de Graauw et al.ll2010l ) 
observations of the ortho- and para-B^O and H^O to- 
wards six of the eight sources in the Guaranteed Time 
Key Program PRISMAS (P.I. Maryvonne Gerin). The 
sources are DR21(OH), G34.3+0.1, W28(A), W33(A), 
W49(N) and W51. Their coordinates, distances and ve- 
locities are listed in Table [T] Figure [T] shows their ap- 
proximate positions in the Milky Way. The observed 
water lines are the three ground state transitions of H2O 
and H^O (2 i2 — l i, lio — loi, and In — 0oo) as well as 
two higher excitation transitions of P-H2O (2o2 — In and 
2n — 202) and one higher excitation transition of 0-H2O 
(2 2 i — 212). Their frequencies, quantum and pseudo- 
quantum numbers, and energies are listed in Table [2j 
The ortho ground state is 34 K above the para ground 
state that we use as a reference for all the energy lev- 
els. Figure [2] shows the water energy diagram with the 
transitions observed with HIFI. 

The data were obtained between March 2010 and April 
2011 using the dual beam switch observing mode. We 
use the HIFI wide band spectrometer (WBS) which pro- 
vides a spectral resolution of 1.1 MHz (0.6 km s _1 at 
557 GHz and 0.2 km s" 1 at 1669 GHz) over a 2.4 to 
4 GHz IF bandwidth. Each H 2 and H^ 8 line has been 
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Figure 1. Position of the sources in the Milky Way. The red and 
blue curves indicate the LSR velocities. The dotted circles show 
heliocentric radii every 3 kpc. The dashed line indicates a Galactic 
longitude of ±90° . 



observed with three different local oscillator frequencies 
(LOF) and two different polarisations (H and V) . All the 
observations have been reduced using hipe (iottii2ninn 
pipeline version 6. For a given H2O or H^O line, three 
different LOFs are not enough to deconvolve and isolate 
a single sideband (SSB) spectrum. Therefore, the H 2 
and H;j, 8 lines are contaminated by features from the 
other sideband. In some cases, these contaminating lines 
fall near the H2O and H^O lines. We do not combine 
the two polarisations immediately, because they target 
positions a few arcseconds away from each other, and 
the receivers and detectors are different, so this may in- 
troduce systematic differences. 

Hereafter, we use W51 as an example. Figure [3] shows 
the double sideband spectra of the nine water lines to- 
wards this source. In the Appendix, Figures [1])] to [30] 
show the double sideband spectra of the nine water lines 
for the five other PRISMAS sources. All the spectra pre- 
sented in this paper show similar properties. Differences 
of about 1 K between the two polarisations are detected 
in the continuum intensity for a few spectra (e.g. W51 
at 752 GHz, W49N at 547 GHz, G34.3+0.1 at 987 GHz). 
We first summarize the detection of absorption and emis- 
sion features and characterize the interstellar clouds. We 
then discuss each transition from the most to the least 
excited in the water diagram. 

2.1. Absorption and emission features 

The six PRISMAS sources we present show H2O 
and H^ 8 in emission at 752, 987, 556, and 547 GHz 
near the systemic velocities of the background sources. 
Some of them also show H2O in emission at 1661 GHz 
(DR21(OH), W28A, W33A), 1669 GHz (DR21(OH), 
W28A, W33A) and 1113 GHz (DR21(OH), W28A, 
W33A, W49N). H^O is seen in emission at 1655 GHz 
(W28A) and 1101 GHz (W28A, W33A). Absorption in 
the ground state transitions of H 2 and H^O is de- 
tected towards all sources though it is sometimes faint 
or difficult to discern at 547 GHz (e.g. towards W51). 
Absorption in the excited transition of H2O is detected 
towards all sources at 1661 GHz and possibly detected 
at 987 GHz. At 752 GHz, no absorption is apparently 
detected towards any of the PRISMAS sources. 

The velocity of the emission components relates to the 
velocity of the source in the Galaxy while that of an ab- 
sorption component depends on the velocity of the line 




Figure 2. Energy level diagram for ortho and para water. The 
transitions within the HIFI spectral coverage are indicated by solid 
lines, while those outside of this range are indicated by dashed lines. 
The transitions observed in PRISMAS are indicated by heavy lines, 
labeled with the frequency of H^O and HJ^O (in italics) when 
observed. 

of sight cloud in the Galaxy. The molecular hydrogen 
column density in the interstellar clouds have been de- 
rived from observation s of CH and HF towa rds the PRIS- 
MAS lines of sight bv lGodard etaLl (f2012h . assuming a 
column density ratio of 0.4 for HF/CH, and an abun- 
dance ratio ra(HF)/n(H 2 ) = 3.6 x 10" 8 . They also give 
the atomic hydrogen column density obtained by various 
authors with the VLA interferometer. The total hydro- 
gen column density of a given cloud is usually between 
a few 10 20 cm~ 2 and a few 10 21 cm -2 though some are 
lower limits. This corresponds to Ay of a fraction of, 
to a few magnitudes. The fraction of molecular hydro- 
gen (2AT(H 2 )/(AT(H) + 2JV(H 2 ))) is typically 30%, but 
varies significantly (from less than 20% to more than 
50%). The interstellar clouds we detect in H2O absorp- 
tion therefore are translucent clouds. The typical density 
and temperature i n those clouds are abo ut 100 cm -3 and 
50 to 100 K (e.g. iRachford et "all [20021) . We note here 
that the water distribution in velo city space has been 
shown to follow that of HF (see e.g. iNeufeld et al.ll2010t 
iSonnentrucker et al.ll2010l ). implying that water, unlike 
OH+ and H2O" 1 " probably arises in the gas that is mainly 
molecular. 

2.2. Excited p-H-z O transition at 752 GHz 

The only water line that does not show any indica- 
tion of absorption is that of P-H2O at 752 GHz. Indeed, 
it is the highest excitation line of water observed with 
HIFI in the PRISMAS progranQ, with a lower level at 
101 K above the para ground state. Therefore, the water 
molecules in clouds along the line of sight, which are be- 
lieved to be relatively unexcited, do not absorb photons 

7 Higher energy transitions are observed with PACS, but not 
velocity resolved. 
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Table 1 

Sources observed. 



Source 


R.A. 


Decl. 


I 


b 


Distance 


Emission 


Absorption 














velocity 


velocity range 




(J2000) 


(J2000) 


(°) 


(°) 


(kpc) 


(km s -1 ) 


(km s" 1 ) 


DR21(OH) 


20 h 39 m 01.P 


+42°19'43" 


81.7 


+0.54 


1.5 


-4 


-5 to 15 


G34.3+0.1 


18 fc 53 m 18.7 s 


+01°14'58" 


34.3 


+0.1 


3.8 


58 


10 to 70 


W28A 


18' 1 00 m 30.4 a 


-24°04'00" 


5.9 


-0.39 


1.9 


10 


-10 to 30 


W33A 


18 k 14 m 39 .4 s 


-17°52'00" 


12.9 


-0.26 


3.8 


37 


20 to 45 


W49N 


19 h 10 m 13.2 a 


+09°06'12" 


43.2 


+0.01 


11.4 


2 


to 70 


W51 


19' 1 23 m 43.9 s 


+14°30'30" 


49.5 


-0.38 


5.5 


55 


to 75 



Note. — We give the equatorial and Galactic coordinates, and distances. We also give the ve l ocity of the water emis sion peak and the 
velocity range over which water absor ption is det ected. References for distances arc Rvgl et al. (2012) for DR21(OH), Fish ct al. I|2003l ) 
for G34.3+0.1. IVeIazquez eTaLl (2002) for W28A, Faund ez~et~aTI pOOl ) for W33A, [Gwinn ct al. (1992) for W49N, and Rcid ct al. j2009ft 
for W51. 



Table 2 

Water lines observed. 





Frequency 




Jk_ 1 k 1 (1) 


Eu/k 


Ei Ik 


gi/gu 


log(A u i) 


Av 


RMS 




(GHz) 






(K) 


(K) 






(m/s) 


(mK) 


(%) 


o-H 2 


1669.905 


2l2 


loi 


114 (80) 


34 (0) 


3/5 


-1.252 


197 


69-357 


2.4-6.5 




556.936 


lit) 


loi 


61 (27) 


34 (0) 


i 


-2.461 


592 


5-44 


0.9-3.6 




1661.008 


221 


2l2 


194 (160) 


114 (80) 


l 


-1.514 


199 


74-348 


2.4-6.8 


p-H 2 


1113.343 


111 


0oo 


53 





1/3 


-1.734 


296 


25-105 


1.6-3.2 




987.927 


2(12 


111 


101 


53 


3/5 


-2.333 


334 


33-133 


2.0-4.2 




752.033 


2n 


202 


137 


101 


1 


-2.152 


439 


19-75 


1.5-4.6 


o-Hf O 


1655.868 


2l2 


loi 


114 (80) 


34 (0) 


3/5 


-1.263 


199 


65-248 


2.4-4.8 




547.676 


lio 


loi 


61 (27) 


34 (0) 


1 


-2.483 


602 


6-27 


1.3-2.8 


p-HfO 


1101.698 


111 


0oo 


53 





1/3 


-1.748 


299 


32-137 


2.2-4.1 



Note. — For each transition we give its frequency, the upper and lower level designations (Jk_ 1 k 1 ), and energies expressed as equivalent 
temperatures relative to the ground state of the para form of that isotopologue (T) , the degeneracy factor ratio (gi/g u ), and the spontaneous 
emission rate (A u i). We also give the ortho energy level relative to their ground state in parentheses and the equi valent velocity reso lution 
for the HIFI spectral resolution of 1.1 MHz in the last column. The spectroscopic data are from the JPL database dPickett et all 19981) . The 
RMS columns give the typical range of the estimated rms in the observed spectra, away from the main absorption and emission features, 
in mK, and relative to the single sideband continuum intensity, in %. 



at this frequency. The asymetry of the emission profile 
indicates the presence of multiple emission components. 
The position of the emission peak defines the velocity of 
the source. We report this velocity in Table Q] for each 
source. These velocities are in very good agreement with 
values from literature. 

2.3. Excited P-H2 transition at 987 GHz 

The next highest excitation line of p-H^O is at 
987 GHz. Its spectrum resembles that of P-H2O at 
752 GHz except that it shows some sign of absorption. 
The absorption feature at 987 GHz is more significant 
towards DR21(OH) than towards G34.3 and W49N, and 
it is weak towards W33A and W51. However, given the 
equivalent temperature of the lower level of this transi- 
tion (53 K), and the apparent agreement in velocity with 
the emission component, we assume this absorption com- 
ponent is due to warm water vapor in the background 
source's envelope or outflow. 

2.4. Excited 0-H2O transition at 1661 GHz 

On the ortho side of the water energy diagram, the first 
excited transition is that at 1661 GHz. Its lower level is 
at a lower equivalent temperature than that of the high- 
est excited transition of the P-H2O at 752 GHz when 



measured relative to the ground state of 0-H2O (80 K vs 
101 K). The 0-H2O transition is always detected in ab- 
sorption and towards three sources in emission. Transi- 
tions of CH at ~1661 GHz, about 100 MHz or 20 km s" 1 
from that of H2O are likely detected towards DR21(OH), 
G34.3+0.1, W49N and W51. We identify the o-H 2 line 
thanks to the matching of its velocity with that of com- 
ponents seen in the spectrum of other H2O and H3, 8 
transitions. Similar to what we observe at 987 GHz, we 
always find a good agreement between the velocity of the 
0-H2O absorption feature and that of the emission fea- 
ture at 752 GHz, which indicates this component arises 
from the background source's envelope or outflow rather 
than from the ISM clouds along the line of sight. 

2.5. Ground state H2O transitions at 1113, 1669 and 
556 GHz 

The three H2O ground state lines, the P-H2O line at 
1113 GHz and the o-H 2 lines at 1669 and 556 GHz, re- 
semble each other to a significant degree for each source. 
They all exhibit multiple absorption features with sim- 
ilar positions, shapes and intensities. Optically thick 
and thin absorption features are found towards each 
source. Some absorption components are wide (e.g. 
~ 20 km s _1 towards W51 at about 60 km s _1 ), some 
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Figure 3. Double sideband water line spectra towards W51. The top row shows the excited transitions of H2O, the middle row shows the 
ground transitions of H2O, and the bottom row shows the ground transitions of H^O. In each panel, we show the H polarisation spectrum 
(black solid line), the interpolated continuum (blue solid line), the best-fit continuum (orange solid line), the emission components (orange 
dashed lines), the best fit (red solid line), the absorption components (red dashed lines), and the zero continuum level for the single sideband 
(grey dashed line). Some of the features seen are not due to water molecules (see text for details). For instance, the 45 km s — 1 emission 
feature at 547 GHz is likely due to 34 S02 and methanol. The 40 km s — 1 absorption feature at 1655 GHz is due to the strongest water 
absorption at 1661 GHz. The 35 and 50 km s _1 absorption features at 1661 GHz are likely due to CH. 



are much narrower (e.g. < 1 km s 1 towards G34 at 
about 27 km s _1 ). We indicate the range of absorption 
velocities for each source in Table [TJ We note here that 
for each PRISMAS source, the velocity and width of the 
emission feature seen at 752 and 987 GHz is such that 
it covers at least partially the range of velocities of the 
absorption components. 

Emission from the wings of water lines is detected at 
556 GHz towards each source, and is visible at 1113 GHz 
and 1669 GHz only towards DR21(OH), W28A and 
W33A. However, given the position, width and intensity 
of the absorption features detected towards G34.3+0.1, 
W49N and W51, we cannot claim there is no water emis- 
sion at those two frequencies. Only a few absorption fea- 
tures are significantly displaced from potential emission 
contamination (e.g. towards G34. 3+0.1 at about 10 and 



30 km s _1 ). All of the absorption features have veloci- 
ties that match those of the Galactic arms crossed by the 
lines of sight, as well as the velocity of the arm in which 
the source is located. However, we cannot claim that all 
features are associated with the foreground ISM and are 
not intrinsic to the source. For instance, the components 
with velocities lower than 30 km s -1 have been suggested 
to be due to W49 itself (see lPlume et alJl2004| ). We de- 
tail the possible association of each absorption feature 
with a Galactic arm below. 

2.6. Ground state H%0 transitions at 1113, 1669 and 

556 GHz 

The three H^O ground state lines (the para transi- 
tion at 1101 GHz and the ortho transitions at 1655 and 
547 GHz) closely resemble each other for a given source. 
The line profiles usually exhibit only one absorption fea- 
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ture at a velocity close to that of the absorption feature 
seen at 987 GHz. An emission feature that resembles 
that observed at 752 and 987 GHz is always detected 
in the spectrum at 547 GHz, though it is significantly 
weaker. Consequently, the spectrum at 547 GHz closely 
resembles that at 987 GHz. This is particularly true 
towards DR21(OH), W28A, and W49N. The emission 
component is also visible at 1101 GHz towards W28A 
and W33A, and at 1655 GHz towards W28A. The abun- 
dance of H^O relative to H2O is generally assumed to 
be about 1/500 though it has been shown to inc rease by 
abou t a factor 2 towards th e inner Galaxy (e.g. iPenziad 
frMll I Wilson fc Roodlll994D . Therefore the H| 8 absorp- 
tion features associated with translucent clouds are not 
detected. 

2.7. Features from other sideband and from other 
molecules 

As mentioned previously, features from the other side- 
band or from other molecules may in a few cases appear 
very close to those of H2O and H3, 8 0. The most common 
case in our analysis is that of the CH (N=2, J=5/2-3/2) 
transition at 1661.1 GHz, at about —20 km s _1 from 
that of H 2 at 1661 GHz. Due to the difficulty of de- 
convolving the two sidebands, a combination of the CH 
and the H2O features also appear in the H5, 8 tran- 
sition at 1655 GHz. Other contaminants include the 
S0 2 28 6 , 2 2-28 5 ,23, 34 S0 2 21 3) i9-20 2 ,i8, and the methanol 
15i,i5l52.i4, Vt = 1 transitions in the H3, 8 spectrum at 
547 GHz' at -70 km s" 1 , +35 km s -1 , and -10 km s -1 , 
respectively, relative to the H^O feature (e.g. towards 
G34.3+0.1, W28A and W49N). 

3. RESULTS 

To estimate the OPR of the water vapor in the translu- 
cent clouds, we need to measure the column density of 
ortho and para water along the line of sight, which in turn 
requires us to determine the optical depth as a function 
of velocity. The main difficulty that we face in estimat- 
ing the optical depth of the water lines is due to the 
background emission. The sources have been chosen for 
their strong dust continuum. However their spectra also 
show emission of water lines. The velocity and width 
of the water emission and absorption components is such 
that they usually overlap significantly, over at least a few 
tens of km s _1 . Therefore, a careful estimate of the wa- 
ter emission lines from the source is required to infer the 
optical depth of the water lines from the clouds along the 
line of sight. This section provides many details on our 
analysis and the reader not interested in those details is 
invited to directly go to the discussion section. 

3.1. Optical depth 

The optical depth is r = — ln(//7o) where I is the 
observed single sideband antenna temperature and Iq 
is the signal sideband pseudo-continuum, which com- 
prises a dust continuum component Idust,sSB and the 
emission features. IdustjSB denotes the dust contin- 
uum component of the image sideband and Idust.DSB 
the dust continuum component of the double sideband 
spectrum. We treat Idust,sSB, I dust , isb and Idust,DSB 
as constants and I dus t,DSB = Idust,ssB + IdustjSB- For 
instance, towards W51, at 556 GHz, Idust,DSB — 2.3 K, 



and Idust.iSB — 1-1 K, as given by the antenna tem- 
perature in the plateau of optically thick absorption fea- 
tures. Since we do not deconvolve the spectra, we treat 
all the emission features as part of the signal sideband 
pseudo-continuum Iq. Towards W51, at 556 GHz, Iq is 
thus represented by all the emission features on top of a 
continuum Idust.ssB — 1.2 K. The actual values of the 
continua depend on the sideband gains gisB and gssB, 
linked by the relation gis b + gss B = !• 

In this paper, we use two methods to measure the 
pseudo-continuum Jo. The first method uses the informa- 
tion provided by the higher excitation transitions (e.g. at 
752, 987 and 1661 GHz) about the water emission compo- 
nents in terms of amplitude, velocity position and width, 
to better constrain the continuum for the lower transi- 
tions. We give details of this method in what follows. 
From this, we obtain a first estimate of the single side- 
band pseudo-continuum intensity io,6/i parameters for 
all absorption and emission components, and the side- 
band gain ratio. The second method consists of interpo- 
lating each spectrum over the velocity ranges where the 
absorption features are detected with a spline function. 
We manually select regions and points to be connected 
for each spectrum. We use the same velocity points for 
both polarisations. Since some spectra may show signif- 
icant noise, we first smooth all the velocity ranges with 
a 10 pixel wide window. With this method, we obtain 
an estimate of the double sideband pseudo-continuum in- 
tensity from which we later subtract the best estimate of 
the image sideband dust continuum determined using the 
sideband gain ratio. We therefore obtain a second esti- 
mate of the single sideband pseudo-continuum intensity 
Ia,int- 

3.1.1. Best fit method 

In this section, we detail the method we use to fit the 
data in order to infer the optical depth of each cloud 
along the PRISMAS lines of sight. For each water tran- 
sition we model the observed intensity as follows: 

I = {gisB 1 1 sb + gssB Issb) , (1) 

gssB 

where: 

gisB=i-gssB, (2) 

IlSB =Idust,ISB, (3) 

Issb = ( Wsss + £ a^-^ 2 ^) 



x exp ( - ^2 a 3 e^ v - v ^' 2 ^ j , (4) 

where the sums over i and j correspond to the emission 
and absorption features, assumed to be gaussians, respec- 
tively. The fitting procedure therefore aims at obtaining 
the best values for gssB, IdustjSB, 1 dust , SSB, and the 
gaussian parameters. 

We start with a fit to the spectrum at 752 GHz as 
only water emission, without absorption, is seen at this 
frequency. This way, we define the properties of the emis- 
sion components (position, width and intensity) . We use 
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as few components as possible - usually two to three - 
to obtain a satisfactory fi t of the spectrum . We use the 
MPFIT package in IDL (|Markwardtll2009t ). We give an 
initial guess specific to each source. 

We then move downwards in energy in the water level 
diagram and fit the spectrum at 987 GHz. We assume the 
properties of the emission components remain roughly 
similar to those defined at 752 GHz: the position of each 
emission component can change by up to 3 km s _1 while 
the width can change by a factor 2. The intensity of 
each emission component is free to vary, as long as it 
remains positive. We then add one or two absorption 
components to fit the spectrum at 987 GHz. For each 
absorption component, we use a single gaussian profile 
to model the corresponding optical depth. We set the 
initial parameters (position, width and intensity) for each 
component manually. As for the emission component, 
we allow the position to be within 3 km s _1 of the initial 
guess and the width to be within a factor 2. 

Then, we move to the H2O ground state transitions, 
where the absorption features due to the clouds along 
the line of sight are detected. We use the spectrum at 
556 GHz to define the absorption components, in addi- 
tion to that already defined at 987 GHz. We again use as 
few gaussian profiles as possible for each of the broad ab- 
sorption components. We usually define one component 
at each local minimum of the observed antenna temper- 
ature. For a given source, we use the results from the fit 
of the spectrum at 556 GHz as the initial guess for both 
the fits at 1113 and 1669 GHz. Therefore, the number 
of emission components and absorption components are 
the same for the three H2O ground state lines. Addition- 
ally, the initial parameters of the components (position, 
width) are identical for the 1113 and 1669 GHz spectra. 
We manually adjust the intensity of the emission com- 
ponents for each spectrum independently. Our method 
is thus consistent with the observed similarity between 
those three spectra. In a few cases at 1669 GHz, we have 
to use additional absorption features to account for con- 
tamination due to other species, either from the same or 
the other sideband (see section |2~7| . 

We finally fit the H^ 8 lines at 547, 1101 and 1655 GHz 
as well as the H2O line at 1661 GHz using the results from 
the fit of the spectrum at 987 GHz as the initial guess. 
As discussed above, the spectra at 547 and 1661 GHz, 
and to a lesser extent at 1655 and 1101 GHz, are very 
similar to that at 987 GHz. We manually adapt the emis- 
sion components intensity to account for the abundance 
difference between H2O and H^O. In a few cases for the 
spectrum at 1655 and 1661 GHz, we have to use addi- 
tional absorption features to account for contamination 
due to other species, either from the same or the other 
sideband (see section |2~T]) . 

We perform all these fits separately for each polarisa- 
tion. However, we use the same initial guess for both 
polarisations observed in a given source and a given wa- 
ter line. We combine the results from both polarisations 
only when deriving the water column density (see section 
1531 . 

3.1.2. Gain and continuum ratio 
8 http://purl.com/net/mpfit 



Since we work on the double sideband spectrum, we 
have to take into account possible variations in the con- 
tinuum due to different sideband gains (gssB 7^ 0.5) 
and to the intrinsic spectral shape of the dust continuum 
{Issb f Iisb). At 752 and 987 GHz, and for the H^ 8 
transitions, this is not critical since there is no optically 
thick water absorption. Therefore, we assume the side- 
band gains to be both equal to 0.5. We have also no way 
to constrain the dust continuum intensities IdustjSB and 
Idust,sSB at these frequencies and in return, they do not 
affect the results of the fit. We therefore let them vary 
within a default range of ±5% away from Idust,DSB /2. 
On the other hand, for the ground state transitions, the 
water absorption features that are optically thick allow 
us to constrain the sideband gain and the dust contin- 
uum. 

We add gssB as a free parameter in the fit for the 
ground state transitions and allow it to vary within 5% 
of 0.5. We justify here these limits. At a velocity where 
the foreground absorption is optically thick, Issb = 
which leads to (1 - gssB) IdustjSB = 9SSB Idust.DSB- 
The limits on the gain value then give 0.9 IdustjSB < 
Idust,DSB < 1.1 IdustjSB- We then compare Idust.DSB, 
as given by the median of the spectrum away from emis- 
sion or absorption features, to IdustjSB, as given by 
the single minimum value of the spectrum where the 
water absorption is optically thick, taking into account 
the standard deviation of the spectrum a. For all our 
spectra, we find that the constraints on Idust.DSB and 
IdustjSB are respected within 3-cr at most. This is very 
satisfying since we estimate IdustjSB with the single min- 
imum value point in each spectrum rather than with an 
average of the intensity in the plateau of optically thick 
absorption features. The values we find for gssB when 
performing the fits are very close to 0.5 and slightly de- 
crease with increasing frequency. When averaging over 
all the PRISMAS sources and both polarisations, we find 
gssB = 0.51 ± 0.02, 0.51 ± 0.01, and 0.50 ± 0.02 at 556, 
1113, and 1669 GHz, respectively. We do not find any 
difference between the two polarisations or any trend 
when comparing the sources, except that gssB is sys- 
tematically the lowest towards W28A with values below 
0.5, whic h may indicate some con tamination by the OFF 
position. Roelfsema ct al. (2012) report in-orbit perfor- 
mance of the HIFI instrument. The measured sideband 
gain ratio is usually between 0.45 and 0.55, through- 
out the entire frequency range of HIFI. Our findings are 
therefore in agreement with expectations. 

We also add the intensities IdustjSB and Idust,sSB as 
free parameters in the fit of the ground state transitions 
to estimate the variation of the dust continuum between 
the two sidebands. We estimate the variations of the dust 
continuum over the entire frequency range covered by our 
data, from 547 to 1669 GHz. We find that the antenna 
temperature follows a law in v 11 , except towards W28A 
where it is closer to v 2A . If we limit this analysis to the 
frequencies below 1600 GHz, the antenna temperature 
follows a law in v 2 3 except towards W28A where it varies 
in f 2 8 . Therefore, at 556 GHz, where the effect is the 
strongest, the difference between the dust continuum in 
the upper sideband and in the lower sideband, which are 
separated by 12 GHz, should be, at the most, about 6% 
((556 + 12) 2 8 /556 2 - 8 ~ 1.06). 

At 556 GHz, the lower sideband contains the water 
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Table 3 

Best-fit gaussian parameters for the absorption features towards 
W51 and column density in the lower level of each transition. 



Transition 




FWHM 


' 


1 y 1 ^ 1 


(GHz) 


(km s ^) 


(km s _1 ) 




(xlO 12 cm -2 ) 


Q87 0H7 


58.1 


7.02 


0.350 


21 ± 2 




67.8 


8.87 


0.264 


13 ± 8 


1661 OH7 


58.6 


6.15 


0.810 


69 ± 14 




66.0 


13.6 


0.442 


77 ± 11 


OOV VjXXlZr 


^ 0*3 


J- .ou 




3 q 4. n c; 
o.y in u.o 




U.O-L 


1 07 




2 i 2 




6.73 


3.18 


713 


10 i 4 




12.6 


3.00 


132 


2 ± 0.5 




45.7 


1.18 


1.84 


10.1 i 0.2 




48.5 


2.51 


0.373 


4 ± 2 




56.9 


5.32 


1.75 


43 ± 2 




65.7 


7.58 


4.97 


185 ± 24 


lllO VTllZ 


^ 08 
o.uo 


1 22 


4^7 






u.^u 


70fi 


^OQ 


8A -t- 0^ 
u.o^t in u.uo 




6.49 


2.62 


0.616 


3.4 it 0.9 




12 8 




04Q7 


Q 4- 1 
u.o in u.j- 






QR1 
u. yu J. 


2 04 


zl fi -t- O 1 




48.5 


2.65 


0.360 


2.2 it 0.2 




56.1 


4.19 


2.44 


24 it 3 




64.9 


6.18 


16.9 


244 ± 25 


1669 GHz 


5.24 


0.896 


0.497 


1.9 it 0.7 




6.58 


2.95 


1.05 


15 ± 1 




6.65 


0.657 


1.04 


3 it 1 




12.4 


3.16 


0.0991 


1.5 it 0.2 




45.8 


1.07 


2.43 


12.1 ± 0.3 




48.8 


3.17 


0.649 


9.6 ± 0.1 




56.4 


4.97 


3.04 


71 ± 4 




65.0 


6.52 


17.6 


534 ± 3 


547 GHz 


57.9 


6.04 


0.157 


4.7 ± 0.9 




64.9 


16.4 


0.141 


10.8 it 0.1 


1101 GHz 


58.1 


8.47 


0.194 


3.7 ± 0.9 




67.6 


4.19 


0.0659 


0.7 it 0.1 


1655 GHz 


58.0 


5.70 


0.318 


9 ± 4 




66.0 


9.23 


0.240 


10 ± 3 



Note. — The uncertainty in the column density in the lower 
level of the transition takes into account the uncertainty in the 
continuum and the standard deviation of the antenna temperature. 
The column densities are derived assuming T ex -C hv/k. 



We here present the results for the sightline towards 
W51. The results for the other lines of sight can be 
found in the appendix. The continuum estimates are all 
shown in Figure 02 All the observations are very well 
reproduced with the best-fit method. The reduced x 2 is 
usually very close to 1. 

The two methods we use to estimate the continuum 
Iq provide us with two ways t o infer the optical depth 
along each line of sight. Figure [4(a) compares the opti- 
cal depth inferred from each method for the three ground 
state H2O transitions. On the one hand, the best- fit 
method provides us with the gaussian parameters for 
each absorption features. We give these parameters in 
Table [3] Uncertainties on the gaussian parameters are, 
in 80% of the cases, lower than 0.2 km s" 1 on the po- 
sition, lower than 0.4 km s _1 on the width, and lower 
than 0.1 on To. The sum of all the absorption gaus- 
sian components gives the optical depth as a function 
of velocity. By construction, the estimate of the opti- 
cal depth with the best-fit method is thus always zero 
at 752 GHz. Also, there is no noise in the inferred op- 
tical depth curves. On the other hand, we combine the 
interpolated continuum with the observed antenna tem- 
perature to compute r = — ln(I/7o,int). Since io,m* is a 
smoothed interpolation of the spectrum, the inferred es- 
timate of the optical depth shows the same noise as the 
data. Consequently, the optical depth at 752 GHz de- 
rived from the interpolated continuum is not zero. Also, 
towards W51, the noise in the data leads to a residual of 
about a few 0.01 in the optical depth derived with the 
interpolated continuum for the ground state transitions. 
We use these residuals to estimate upper limits for the 
column densities of undetected absorption features (e.g. 
excited transitions of H2O in translucent clouds). 



3.2. Column density 

We compute the column density Ni of H2O and H3, 8 
in a given lower energy level using: 



features therefore we expect Idustjss > Idust,ssB- We 
thus assume IdustjSB is about 5% larger than the mean 
value of the continuum, as given by Idust.DSB /2, and 
Idust,ssB is about 5% smaller. We allow both param- 
eters to vary by up to 5% around those values, which 
means that in the extreme cases, they can be either equal 
or about 20% apart from each other. At 1669 GHz, the 
upper sideband contains the water features therefore we 
expect Idust,sSB > IdustjSB- We thus use the same 
limits as defined above but on opposite parameters. At 
1113 GHz, we combine data from lower and upper side- 
bands therefore we expect IdustjSB — Idust,ssB- We 
thus let vary the parameters within ±5% of Idust,DSB /2. 
The results from the fit confirm our expectations with 
hustjSB /Idust,ssB = 1.05 ± 0.05, 1.01 ± 0.05, and 
0.91 ± 0.05 at 556, 1113, and 1669 GHz, respectively. 
As for gssB we do not find any systematic between po- 
larisations and we find that apart from W28A, for which 
I dust , isb I Idust,ssB is almost systematically the highest, 
there arc no obvious trends when comparing sources. 

3.1.3. Results 



N, 



,9i v 



-Av 



g u c 3 A ul i_exp(^) 



(5) 



where gi/g u is the ratio of the degeneracy factors of the 
lower and upper levels, v is the water line rest frequency, 
A u i is the spontaneous emission rate of the transition, 
T ex is the excitation temperature, and Av is the width 
of the feature. The transition parameters are given in 
Table [5J If we assume T ex hvjk, then the column 
density is given by: 



Nt 



g u c A u i 



(G) 



in SI units. Therefore the corresponding expressions for 
the column densities of each H2O transitions, in units of 
10 12 cm -2 , with Av the FWHM of a gaussian profile in 
km s , and the associated conditions on the excitation 
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temperature are: 

iV 0oo = 2.33 t(1113 GHz)Av for T ex < 53 K, (7) 

N lol = 4.66 r(556 GHz)Av for T ex < 27 K, (8) 

N lm = 4.66 r(1669 GHz)Av for T ex < 80 K, (9) 

AT Xii = 9.23 r(987 GHz)Av for T ex < 47 K, (10) 

iV 2o2 = 5.63 r(752 GHz)Av for T ex < 36 K, (11) 

iV 2i2 = 13.9 r(1661 GHz)Av for T ex < 79 K, (12) 

Similarly, for the H3, 8 transitions, we have: 

iV 0oo = 2.33 r(1101 GHz)Av for T ex < 53 K, (13) 
N lol = 4.66 r(547 GHz)Aw for T ex < 26 K, (14) 
AT lol = 4.66 r(1655 GHz)Av for T ex < 79 K, (15) 

Since we have two estimates of the optical depth, de- 
rived from two estimates of the continuum, we infer two 
estimates of the column density. In Table [3J we give the 
water column density in the lower level of each transi- 
tion for each gaussian absorption component using the 
gaussian parameters (tq for t and FWHM for Av) from 
the best fit. The column densities that we report are 
the average of the two polarisations. The absorption 
features associated with the source (i.e. outflow, enve- 
lope) are always optically thick for the ground state H2O 
transitions and are the only ones detected in the excited 
transitions of H2O and the ground state transitions of 
H-^O. The optical depth and H2O column density de- 
rived for the ground rotational levels of these absorp- 
tion features are thus unreliable while those of the ex- 
cited transitions of H2O and the ground state transitions 
of H^O are reasonable, though usually highly uncer- 
tain. On the other hand, the absorption features associ- 
ated with the translucent clouds are only detected in the 
ground state H2O transitions and not in the H2O excited 
transitions or ground state H^O transitions. Therefore, 
for these absorption components, by construction, the 
column density in the H2O excit ed lev els and the H3, 8 



ground state levels is zero. Figure 4(b) shows the column 
density curves for the ground states transitions towards 
W51 for a common velocity resolution of 0.5 km s -1 . 
In Table 2] we then give the total water column density 
iV(H 2 0) = iVooo + N Un + N Ul + N 2o2 + N 2l2 for each 
absorption component identified by the velocity position 
of its peak. We only compute the total column density 
of water for the absorption components associated with 
translucent clouds. The total column density of water 
for the absorption component associated with the source 
is not only highly uncertain but also not representative 
of the true column density of water, which also com- 
prises the molecules that are responsible for the emission 
features and that we do not constrain. For the interstel- 
lar absorption features, Ni lt , N 2o2 , and N 2l2 are zero 
by construction. Therefore -ZV(H 2 0) = iVi 01 + Nq 00 and 
the fractional abundance of H2O in the ground states, 
/ S (H 2 0) = (N lol + N 0oo )/N(U 2 0), is thus equal to 1. 
The column density in the ortho ground state iVi 01 is 
given by two H2O transitions (556 and 1669 GHz for 
H 2 0). We show in the discussion (see section |4"T21 that 
the 1669 GHz transition gives a significantly better esti- 
mate of Ni 01 than the 556 GHz transition. The values 
presented in Table [4] and in the Appendix take this into 



account. In most cases, the total H2O column density is 
about a few times 10 12 to a few times 10 13 cm~ 2 . 

In order to estimate more accurately the possible con- 
tribution of the excited transition to the total column 
density of water, we convert the optical depth curves de- 
rived from the interpolated continuum into column den- 
sity curves. Since a velocity bin is not a gaussian profile, 
the column density inferred from Equation [6] needs to be 
multiplied by a correction factor 2-^/2 ln(2)/v / 27r ~ 0.94 
and Av is given by the velocity bin width. Figure 
|4(b)| compares the inferred column density curves for the 
ground states transitions towards W51 for a common ve- 
locity resolution of 0.5 km s _1 to that previously derived 
from the gaussian components. For the ground state H2O 
transitions, the two estimates of the column density are 
very similar at the position of the absorption features. 
Away from these absorption features and for the other 
transitions (excited transitions of H2O, ground state 
transitions of H^O), the column density curves only 
show the residual due to the small difference between 
/ and lo.mt- With these curves, we compute the water 
column density of each lower level, for each absorption 
feature, by integrating the column density curves over the 
velocity range that best covers the absorption feature. A 
given absorption feature usually combines several gaus- 
sian components, as shown in Table |H For the excited 
levels of H 2 (i.e. 2Vi 115 N 2o2 , and N 2l2 ) and the ground 
states of H^O the derived values are upper limits. They 
usually are about a few 10 11 cm -2 . For the ortho ground 
state level of H2O (i.e. I01), we use the column density 
inferred from the transition at 1669 GHz rather than 
that at 556 GHz, for the reasons detailed in the dis- 
cussion (see section |4~2|) . In Table IH we give the total 
water column density over several velocity ranges. For 
the translucent clouds, the total column density N(H 2 0) 
is thus given as iVo 00 + 2Vj. 01 with a lower uncertainty 
S(Nq 00 +-ZVi 01 ) and an upper uncertainty 6(Nq 00 +Ni al ) + 
Ni 1± +N 2q2 +N 2l2 that shall be viewed as an upper-limit. 
For each velocity range, we then compute the lower-limit 
of the fractional abundance of water in the ground states 
/ 9 (H 2 0) = min(A l01 + i\T 0oo )/max(iV(H 2 0)). We use 
min(A lol + N 0oo ) = A 0oo + N lol - d(N 0oo + N loi ) and 
max(JV(H 2 0)) = min(A^ lol + A 0oo ) + N Ul + N 2(l2 + A 2l2 
as we have to use a consistent value for the column den- 
sity in the ground states. We also give the upper-limit of 
the H3, 8 total column density by adding the upper lim- 
its of 7Vo 00 and Ai 01 , where AT lol is the average from the 
547 and 1655 GHz spectra, weighted by their respective 
uncertainty. The column density curve also allows us to 
compute the water abundance relative to that of H 2 in- 
ferred from CH and HF absorption lines bv lGodard et al.l 
(2012). The velocity ranges we use to compute the water 
column density are usually very close to those obtained 
from the CH and HF obser vations. We therefore give the 
column density of H2 from iGodard et al.l (|2012[ ) for the 
closest velocity range from their Table 4. We chose to 
report the values they derive from the CH observations 
(their column d) as they saturate significantly less often 
than those derived from the HF observations. We finally 
compute the water abundance for each velocity range. 
We discuss the water column density and abundance in 
section O to O 
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Figure 4. Optical depth and column density in the lower level for the ground state H2O transitions towards W51. In each panel, the 
black solid line is inferred from the gaussian components and the dashed colored line is derived from the interpolated continuum. The 
optical depth is shown at the native velocity resolution of the spectrum while the column density is shown with a common velocity bin of 
0.5 km s — . The top, middle and bottom panels show the 556, 1113, and 1669 GHz transitions, respectively. 



Table 4 

Column density, OPR, and abundance of water towards W51 for the translucent clouds. 





Gaussian decomposition 








Column density curves 






V 

(km s 


iV(H 2 0) 
- 1 ) (10 12 cm" 2 ) 


OPR 


V 

(km s" 1 ) 


AT(H 2 0) 
(10 12 cm" 2 ) 


/ S (H 2 0) 


OPR+ OPR 9 


AT(H 2 ) 
(10 20 cm" 2 ) 


X(H 2 0) 

(10- 8 ) 


(10 12 cm~ 2 ) 


5 
6 
7 


2.8 ± 0.9 
3 ± 1 
18 ± 2 


2 ± 2 

3 ± 1 

4 ± 2 


3 to 15 


27 3+ 4 


>0.88 


3.0±0.8 3.5±0.2 


4.3±1.4 




<0.9 


13 


1.9 ± 0.2 


4.6 ± 0.7 
















46 
49 


16.7 ± 0.3 

11.8 ± 0.2 


2.7 ± 0.1 
4.3 ± 0.5 


44 to 51 


30±? 


>0.49 


5±2 3.2±0.1 


2.6±0.8 


12±f 


<0.8 



Note. — The left hand side of the table gives values derived from the gaussian decomposition of the optical depth while the right hand 
side gives values derived from the optical depth curve. The upper hand error bars on iV(H20) and X(H 2 0) are derived from upper limits 
on the column density in the excited levels and should also be viewed as upper limits. The molecular hydrogen column densities ./V(H2) 
are derived from CH observations ass uming a column density ratio of 0.4 for HF/CH, and an abundance ratio n(HF)/n(H 2 ) = 3.6 X 10" 8 
IjGodard et al. 20l2). See section [3.31 for details about the OPR values and the meaning of the error bars. 



3.3. Ortho-to-para ratio 

We compute the ortho-to-para ratio OPR = (N\ 01 + 
N 2l2 ) / (N Qoo + N lll + N2 02 ). Using the gaussian decom- 
position of the optical depth, for translucent clouds, the 
OPR simplifies into OPR — Ni 01 /N 0oo . We report this 
value as OPR 9 in Table |4] Using the column density 
curves, we compute two extreme values of OPR assum- 
ing (1) there are no water molecules in the lower energy 
levels of the transitions not observed (i.e. for instance 
Ni 10 , -^2 2 u an d N211 are null) and either (2a) the col- 
umn density in the ortho excited levels is null while that 
in the para excited levels is exactly the upper limit (i.e. 
OPR = N lm /(N 0oo + N Ul + N 2o2 )), or (2b) the column 
density in the para excited levels is null while that in 
the ortho excited levels is exactly the upper limit (i.e. 
OPR = (N loi + N 2l2 )/N 0m ). We define OPR+ as the 
mean of these two values. Half their difference is added 
to the uncertainty on OPR 9 to define the uncertainty on 
OPR + . We report this very conservative value in Table 
[31 By defintition, OPR 9 is thus included in the range 
defined by OPR+ ± dOPR+. 



3.4. Uncertainties 



In the following, we detail our estimation of the uncer- 
tainties in the optical depth, column density and OPR. 
Every time we average two estimates of a value (e.g. col- 
umn density from two polarisations), we take into ac- 
count the respective uncertainty into account, using a 
weighting equal to the inverse of the respective uncer- 
tainty squared: 



<X>=5>,A7f)/]T(l/a, 2 ). 

i i 

The uncertainty in the weighted mean is: 
£?((*<- <X» a /°?) 



a 2 = 



in- 1) e:(im 2 ) 



(16) 



(17) 



Every time we add values (e.g. to compute the total wa- 
ter column density), we add the uncertainties in quadra- 
ture: a 2 = £>f). 

One source of uncertainty in r is due to the noise SI 
in the observed antenna temperature /, while the second 
source of uncertainty is the difference <5io between the 
two estimates of the continuum / ,b/ an d Io,int- <^o is 
not constant for a given spectrum and usually is about 
a fraction of Kelvin. Therefore, SIq/Iq usually is about 
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a few percent. To estimate 51, we measure the stan- 
dard deviation of I away from emission and absorption 
features. Consequently, 51 is constant for a given spec- 
trum. At 547 and 556 GHz, 51 ranges from 0.006 K to- 
wards DR21(OH) and W33A, to 0.027 K towards W51. 
At 1655, 1661 and 1669 GHz, 51 ranges from 0.07 K to- 
wards W33A to 0.350 K towards G34.1+0.1. There is a 
slight difference, usually less than a few percent but up 
to ^25%, between the value of 51 inferred from each po- 
larisation spectrum, though we do not see any systematic 
effect. The characteristic value for 51 /I is about a few 
percent. 51 and 51 /I can later be significantly reduced by 
working at a common lower velocity resolution. We chose 
to work with a velocity bin size of 0.5 km s _1 which is 
slightly smaller than the worst resolution of ~ 0.6 km s _1 
at 556 and 547 GHz. 

We perform the best fit to the observations taking the 
uncertainty 51 into account. The resulting parameters 
for the gaussian components (to, FWHM and vq) and 
their uncertainty are computed by the MPFIT pack- 
age. The uncertainty in the column density then is 
5N 0oo = 2.33 x (5t x FWHM + 5FWHM x r ) and 
similarly for each transition and each polarisation sepa- 
rately. We then add the contribution of 5Io. The cor- 
responding uncertainty in the optical depth is 5r con t = 
( T bf ~ r int) — In {Io,bf/Io, int) which is not constant as a 
function of velocity. Therefore, for each gaussian compo- 
nent, we compute the mean value of 5T con t over a veloc- 
ity range centered on and as wide as the FWHM. We 
then convert 5T cont into 5N 0oo = 2.33 x 5T cont x FWHM 
and similarly for each transition and each polarisation 
separately. We then quadratically sum both sources of 
uncertainty in N and average both polarisations. We av- 
erage the gaussian parameters the same way. We report 
the resulting values in Table [3] In Table |4j we give the 
total column density and average OPR of water for each 
gaussian component. 

In Table 21 we give the results from the integration of 
the column density curves, then used to compute the 
OPR and water abundance relative to H2. The un- 
certainty associated with the column density curve is 
5N = N/t x (51 /I + 6I0/I0) for each transition and each 
polarisation separately. It is therefore not constant as a 
function of velocity. When integrating over a given ve- 
locity range, we quadratically sum the uncertainty asso- 
ciated with each velocity bin before averaging the results 
from the two polarisations. When averaging the column 
density derived from each polarisation, we note that in 
the very large majority of cases, the difference between 
the two values is lower than the uncertainty and almost 
always less than 3 times the uncertainty. 

3.5. Comparing the two methods 

We measure the total column density of water us- 
ing the gaussian profiles and using the column density 
curves. Differences between the two methods may arise 
from contamination due to non- water absorption feature, 
possible water in the excited levels, nearby water absorp- 
tion feature associated with the clouds, or broad water 
absorption profiles associated with the source. For in- 
stance, when measuring the total c olumn density towards 
W51 with the curve (see Figure |4(b)| , the feature at 
49 km s _1 is obviously contaminated by the interstellar 
feature at 46 km s _1 and by the broad feature between 



55 and 80 km s" 1 associated with the source. Addition- 
ally, when integrating the curve, we compute the possi- 
ble contribution from excited levels, assuming these are 
upper limits. Finally, the use of finite limits of integra- 
tion may underestimate the column density inferred from 
the curve. Consequently, there might be some scattering 
and/or systematics between the total column density of 
water derived with each method. 

Figure [5] compares the values obtained for each group 
of absorption features with the gaussian profiles and the 
column density curves. In sum, we have 13 groups of 
absorption features, each of which c ould b e an isolated 
cloud or a group of clouds. In figure 5(a) we show that 
over almost two orders of magnitude, the total column 
densities inferred from each method are very well corre- 
lated. The best linear fit to the data leads to a slope of 
1.06±0.01, if we only take into account the uncertainty 
in the ground state population. The slope is 1.03±0.04 
if we add the upper limits on the excited states popu- 
lation to the total column densities. This means there 
is a good agreement between the two methods but that 
a slight systematic effect may lead to lower value of the 
total column density of w ater w hen using the gaussian 
components alone. Figure 5(b) then shows the ratio of 
the column density derived from the two methods as a 
function of the total column density of water. We de- 
tect no clear trend though larger column densities seem 
to give a better agreement between the two methods. A 
linear fit leads to an intercept of ~ 1.1 and a slope of 
- 10~ 3 . The average ratio is 1.05±0.01 or 1.03±0.04 de- 
pending on the uncertainties we use. We conclude that 
there is a systematic effect, whose origin may be one of 
the explanations given above. However, we also show 
here that the differences between the two methods are 
not significant (less than 5%). Both methods thus ap- 
pear to yield equally reliable measurements of the total 
column density of water. 

4. DISCUSSION 

In this section, we discuss in detail the water column 
density, abundance, and OPRs found along the PRIS- 
MAS lines of sight. We first discuss the presence of water 
in the excited levels and quantitatively analyze the dif- 
ferences between the two ortho ground state transitions 
to derive constraints on the physical conditions within 
the translucent clouds. We then summarize our findings 
regarding the total water column density, its abundance 
relative to molecular hydrogen, and discuss possible vari- 
ations throughout the Galaxy. We finally analyse the 
OPR, and how these relate to the physical conditions 
the H2O molecules have interacted with in their history. 

4.1. On the abundance of H%0 and the detection of 
Hi O in the excited transitions 

For the first time, Herschel/HIFI enables the observa- 
tions of all the HJ^O ground state transitions as well as 
the first few excited transitions of both ortho- and para- 
H2O. The only absorption features that are unequivo- 
cally detected in those transitions are related to the back- 
ground source and not discussed further in this paper. 
The absorption features associated with the translucent 
clouds only lead to upper limits of the total column den- 
sity of H3, 8 and to that of the column density of H2O 
in the excited transitions. 
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Figure 5. Comparison between the total column density of H2O inferred from the gaussian profiles and from the column density curves. 
Each data point represents one absorption feature group. In each panel, the red solid line is the best linear fit to the data, the blue dashed 
line is a column density ratio equal to unity, and the orange dashed line is the average column density ratio 



The abundance ratio between 16 O and 18 O is com- 
monly accepted to be ~500, though significant variations 
have been observed through the Galaxy. Therefore, the 
optical depth of H^O is ~500 times lower than that of 
H2O for any given velocity bin. This implies that cither 
the H^O optical depth is too small to be detected or the 
H^O optical depth is large and the absorption feature is 
saturated. Therefore, we can only determine the column 
density of one of the isotopic form of water at a time and 
consequently, can not derive the H^O-to-H^O isotopic 
abundance from our data. For the translucent clouds, 
we only measure iV(H20) and provide upper limits on 
7V(H^ 8 0), usually around 10 12 cm . The inferred lower 
limits on the 16 0-to- 18 abundance ratio are always con- 
sistent with, though significantly lower than 500. 

The excited transitions of water are all detected in 
emission and, with the exception of the 752 GHz para- 
H2O transition, in absorption at the velocity of the 
source. The translucent clouds only exhibit absorption 
features in the ground transitions and only upper limits 
on the column density of water in the excited transitions 
can be computed. Therefore we can estimate the frac- 
tion of the population of water molecules that lies in the 
ground state. We define f g = (N 0oo + Ni 01 )/(N 0oo + 
Ni 11 + N 2o2 + Ni 01 + N 2l2 ) for each absorption compo- 
nent. We therefore assume there are no water molecules 
above the excited levels observed with Herschel. For the 
gaussian components, f g — by construction as no gaus- 
sian are required to fit the excited transition spectrum for 
the translucent clouds. For the column density curves, 
we compute the upper limits of iV2 02 , iVi 01 , and N 2l2 
and use them to derive a lower limit on f g (see Table 
[4j. Towards the six PRISMAS sources, the lower limit 
on f g varies between 0.50 and 0.98. Averaged over all 
the interstellar features where we measure it, the lower 
limit on f g is 0.88, if weighted by the total water column 
density of each absorption feature. We also compute the 
fractional abundance of water in the ground state for the 
ortho and para form of H2O independently for each ab- 
sorption feature and do not find any significant difference 



between them though we notice a better correlation be- 
tween f g and the ortho fractional abundance than the 
para fractional abundance. 

These constraints on the population in the excited lev- 
els thus enable us to estimate the excitation temperature. 
We have: 
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for the para water. We assume here that a single T ex can 
properly describe the water populations we observe in the 
first few ortho or para levels but treat the ortho and para- 
H2O independently. We then subsitutc these expression 
of the column density ratios into the expression of f a and 
f p , the fractional abundances in the ground state levels of 
the ortho and para water, respectively. We thus infer the 
following relations between the excitation temperature 
and the fractional ground abundances: 
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Since the fractional abundances are lower limits, the 
excitation temperatures we compute are upper limits. 
The para-H20 transitions are more sensitive to T ex as 
the para 1113 GHz transition corresponds to a smaller 
temperature difference (53 K) than the ortho 1669 GHz 
transition (80 K). As a consequence, the upper-limit on 
the excitation temperature derived from the para water 
is almost always more stringent. The upper limits on 
the excitation temperature we find for the 13 translu- 
cent clouds range from about 13 to 26 K. The best situa- 
tion corresponds to the cloud at 16 — 28 km s _1 towards 
W28A. If instead of assuming a single T ex for the first 
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three para levels, we compute it using Ni 1± and Nq 00 , we 
find almost exactly the same results, since the correction 
due to the population in the 2q2 level is not significant. 

4.2. On the discrepancy between the two ortho ground 
state transitions 

While SWAS observed the ortho-water ground state 
transitions at 547 and 556 GHz and ISO observed 
the ortho-water ground state transitions at 1655 and 
1669 GHz, the HIFI instrument on board Herschel allows 
us for the first time to observe all these transitions with 
the same instrument. If the assumption of low excitation 
temperature (hv 3> kT ex ) is verified, then our estimates 
of iVi 01 inferred from either the 556 or the 1669 GHz 
transition for H2O should be the same. 



Figure |6(a) shows the column density of the loi level 
inferred from the 556 GHz transition (iVi 01 (556 GHz)) 
as a function of that inferred from the 1669 GHz tran- 
sition (iVi 01 (1669 GHz)) for each gaussian component, 
not associated with the source, in the best fit of the 
spectra. We also show a 1-to-l correlation and the best 
linear fit to the observations. While the two estimates 
of the column density are very well correlated, there is a 
clear systematic effect. The 556 GHz transition statisti- 
cally leads to a lower column density than the 1669 GHz 
transition. The best linear fit has a slope of 0.94 ± 0.01 
and an intercept of —0.2 ± 0.3. We do the same com- 
parison with the values derived from the column den- 
sity curves and find a slightly better correlation (not 
shown here) and N\ Q1 (556 GHz) re mains statistically 
lower than A r i 01 (1669 GHz). Figure 6(b) then shows 
the ratio between the column density interred from each 
transition as a function of iVi 01 (1669 GHz) to look for 
second order variations. The discrepancy is even clearer 
on this plot and seems to increase as the column density 
increases, though the scatter is significant at all column 
densities. A linear fit to the data gives a slightly de- 
creasing slope of ~ 10 -3 and an intercept of ~ 0.9. The 
average value of iVi 01 (556 G-ffz)-to--/Vi 01 (1669 GHz) is 
between 0.87±0.03 and 0.89±0.04 whether we use the 
gaussian components or the column density curves. We 
draw the reader's attention to the fact that the plots 
only show the gaussian absorption component associated 
with the translucent clouds but that some of these might 
be more affected than others by the emission from the 
source. We removed the few most extreme data points 
before running the fit. It is however clear that a system- 
atic discrepancy of about 10% exists between the two 
estimates of the ortho ground state column density. To 
interprete this discrepancy, we explore several possibili- 
ties. 

First, the assumption on the low excitation tempera- 
ture (see equation ([6])) may not be always true. There- 
fore, in the general case, we do not find T556 = T1669 
because r 7^ In {I/Io) but: 



In 



I — B u (T ex ) 
4 — B v (T ex ) 



(21) 



where B v is the Planck function. The condition T ex <C 
hv/k is more easily verified at 1669 (hv/k ~ 80 K) 
than 556 GHz (hv/k ~ 27 K). Additionally, the dust 
continuum is usually significantly higher at 1669 than 
556 GHz. In our data, we observe that the continuum 



intensity is between 6 and 16 times higher at higher fre- 
quency. As a consequence, the optical depth inferred 
from the 1669 GHz spectrum assuming a low T ex should 
be closer to the optical depth derived without any as- 
sumption on T eX7 than that derived from the 556 GHz 
transition. To confirm this and to estimate T ex in the 
translucent clouds, we look for the value of T ex that gives 
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where 4,556 and 4,1669 are the continuum intensities as 
inferred from the best-fit at 556 and 1669 GHz, respec- 
tively, and T556 and T1669 are the optical depth we derived 
under the assumption T ex <C hv/k. In this equation we 
can neglect the contribution of the CMB whose antenna 
temperature at 556 GHz is less than 2 mK. We also as- 
sume here that no other effect can produce any discrep- 
ancy between T556 and Ti669 (see below for other explana- 
tions). For all six sources, over regions where the optical 
depth is larger than 0.1 and smaller than 2.5, we then 
find that the averaged value of T ex ranges from 4.5 K 
towards W33A to 6.7 K towards W51. The required T ex 
varies within each cloud (see Figure [7]). These expected 
values on the excitation temperatures are in agreement 
with the constraints derived from the upper limits on the 
column density of water in the excited levels (see section 
14. lj) though they are more stringent. Assuming that col- 
lisional excitation is responsible for populating the ex- 
cited levels, and that the excitation temperature is 5 K, 
the fractional abundance of water in the excited levels 
(l-/ ) and (1-fp), as derived from Equation[20l would be 
about 2 10~ 7 for ortho-H 2 and about 8 10" 5 for para- 
H2O. These constraint are significantly lower than the 
upper limits derived from the observations of the excited 
transitions. 

The constraints on the excitation temperatures pro- 
vide us with a constraint on the density in the translu- 
cent clouds alon g the line of sight. We us e the computer 
program Radex (jvan der Tak et alJ 120071 ) and the asso- 
ciated Leiden Atom ic and Molecular Database (LAMDA 
ISchoier et al.l 12005) to compute the excitation temper- 
ature of ortho- and para-H ^O. The collision rate s be- 
tween H 2 Q and H 2 are from iDubernet et~aLl (|2009h and 
IDaniel et al.l(l20H . We use a kinetic temperature of the 
gas below 100 K and densities of H2 lower than 10 6 cm~ 3 . 
We assume the only partner for collisions is molecular hy- 
drogen. The results, not shown here, indicate that for a 
kinetic temperature of 50 to 100 K, typical in translucent 
clouds (jRachford et al.l[2002l ) . an excitation temperature 
of 5 K is found for density smaller than ~ 10 4 cm~ 3 for 
ortho-H20. These upper limits on the density are consis- 
tent with, though significantly larger than the densities 
that characterize translucent clouds (~ 100 cm -3 ), and 
do not rule out the presence of high density clumps in the 
clouds. The total hydrogen column densities are about a 
few 10 21 cm -2 while the sizes of the clouds are about a 
few 100 pc at most (see section |4~2|) . which corresponds 
to average densities of a few cm -3 at least. However, the 
cloud sizes are overestimated due to likely velocity dis- 
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Figure 6. Comparison between the column density of the loi level inferred from the 556 GHz transition (iVi 01 (556 GHz)) and the 
1669 GHz transition {Ni 01 (1669 GHz)). Each data point represents one absorption gaussian profile. In each panel, the red solid line is the 
best linear fit to the data, the blue dashed line is a column density ratio equal to unity, and the orange dashed line is the average column 
density ratio. 




-40 -20 20 

Velocity (km/s) 

(a) Towards DR21(OH) 



20 40 60 SO 10' 
Velocity (km/s) 

(b) Towards G34.3+0.1 



50 50 

Velocity (km/s) 

(c) Towards W28A 





20 40 60 
Velocity (km/s) 

(d) Towards W33A 



-20 20 40 6 
Velocity (km/s) 

(e) Towards W49N 




20 40 60 
Velocity (km/s) 

(f) Towards W51 



Figure 7. Comparison between the optical depth derived from the 556 GHz (in black) and the 1669 GHz (in blue) transitions under the 
low T ex assumption. The values of T ex (right y-axis) required to make these optical depths consistent are indicated by the orange crosses 
over the regions where the optical depth is between 0.1 and 2.5. 



persion within them and the inferred density is averaged 
along the line of sight while significant variations may 
occur within each cloud. 

We also compare the corrected n 0I to that derived 
under the assumption of low T ex for both ortho ground 
transitions. The correlation is almost perfect with Ti669 
(see Figure IU). This plot confirms that the 1669 GHz 
transition is better than the 556 GHz transition for esti- 
mating the column density of water in the ortho ground 
state. We discuss how this constraint on the excitation 
temperature translates into a constraint on the density 
in section |4~TI 



Another interpretation of the discrepancy between T556 
and T1669 relies on the significant amount of water emis- 
sion at 556 GHz which makes it more difficult to properly 
constrain the continuum emission and consequently the 
optical depth and column density for this transition. For 
instance, towards W51, the peak of the water emission 
is at least 1 K above the ~1 K continuum at 556 GHz 
and almost undetectable above the ~7 K continuum at 
1669 GHz. The differences in the emission features conse- 
quently affect the spectral shape of the absorption com- 
ponents at 556 and 1669 GHz (e.g. towards W51 be- 
tween 45 and 65 km s _1 ). The spectral shape is also af- 
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Figure 8. Corrected optical depth for the ortho ground state 
of water as a function of the uncorrected optical depth derived 
from the 1669 GHz transition. Each data point represents one 
velocity channel where the corrected optical depth is above 0.1. 
All the sources are combined. The dashed blue line is the 1-to-l 
correlation. 

fected by the velocity resolution (~ 200 m/s at 1669 GHz, 
~ 600 m/s at 556 GHz) and the sensitivity at each fre- 
quency. However, these differences should not introduce 
systematics. The only explanation for the observed sys- 
tematics would be an underestimation of the water emis- 
sion from the source at 556 GHz, or an overestimation 
of the water emission from the source at 1669 GHz. The 
latter seems difficult to obtain as the two estimates of 
the continuum, from the best fit and the interpolation, 
tend to minimize the contribution from water emission. 
The former seems more plausible as the only constraints 
on the water emission profile are usually limited to their 
wings and sometimes a few velocity channels significantly 
above the dust continuum. A perfect example of this 
situation is seen towards W51. The true contribution 
of the water emission to the continuum could obviously 
be much stronger than that derived from the best fit 
and the interpolation at 556 GHz but it could hardly be 
fainter at 1669 GHz (see Figure©. The iVi 01 (556 GHz)- 
to-JVi 01 (1669 GHz) ratio for the two absorption compo- 
nents towards W51 is 0.88 from 3 to 15 km s" 1 and 
0.70 from 44 to 51 km s^ 1 . Therefore, to explain these 
discrepancies with an increase of the water emission at 
556 GHz, the optical depth should be 1.1 and 1.4 times 
larger than that derived from the best fit for each absorp- 
tion component at this frequency. We convert these scal- 
ing factors for the optical depth into correction factors 
for the continuum intensity. From 44 to 51 km s _1 for 
instance, the corrected continuum l§^ corr would be given 
by Io,corr/I = (Io/I) 1A - Averaged over the width of the 
feature, this implies a continuum about 0.3 K brighter. 
This correction seems reasonable though a bit large com- 
pared to the amplitude of the dust continuum (~1 K) 
and that of the water emission inferred from the best- 
fit (~1 K at the peak). It thus seems unlikely that an 
underestimation of the continuum at 556 GHz alone ex- 
plains the discrepancy observed between the two ortho 
ground state transitions. 

One last interpretation we explore is related to the dif- 
ference in the beam size of the observations. The half 
power beam width of HIFI, in band 1 (480-636 GHz) is 
R" in band 6 (1426-1703 GHz). At least 



two effects could be at play. Either the emission from the 
background source or the absorption from translucent 
clouds along the line of sight could vary significantly over 
angular scales comparable to the beam size. Therefore, 
the continuum averaged over the beam could probe dif- 
ferent environments in the source at 556 and 1669 GHz. 
This would in turn affect the flux-to-continuum ratios 
and explain the discrepancy between the two estimates 
of 2Vi 01 . To explain the discrepant values between the 
two ortho ground transitions, one would need a stronger 
continuum at 1669 GHz than at 556 GHz, which can 
be explained by a stronger dilution at the shorter fre- 
quency. To better understand how the beam size af- 
fects the observations, we compare our dat a to those ob- 
tained with SWAS at 5 56 GHz (see e.g. iNeufeld et all 
I2002t iPlume et all 12004 towards W51 and W49 respec- 
tively). The absorption components associated with the 
background source are narrower and shallower in the 
SWAS data than in the HIFI data. In particular, to- 
wards W51, the SWAS observations at 556 GHz show 
only one strong, ~10 km s _1 wide, absorption feature at 
65 km s _1 while our HIFI observations show a very broad 
(~20 km s^ 1 ) feature with two clearly distinct compo- 
nents. Towards W49, the absorption features between 
and 20 km s _1 are well separated in the SWAS data 
while they are blended in the HIFI observations. The 
emission features also seem to change in width and in- 
tensity, though it is more difficult to assess based on their 
wings only. This means that within the SWAS beam of 
3.3' by 4.5', at least some of the absorption features due 
to the outflow or envelope of the source are significantly 
diluted. Therefore, the angular size of these components 
is notably smaller than ~4'. We can therefore assume 
the characteristic scales for the emission components in 
the background are also significantly smaller than ~4' 
which could lead to sensible variations in the continuum 
between the observations at 556 and 1669 GHz. This 
is also suggested by the difference in the continuum's 
intensity between the observations made in the two po- 
larisations, which are usually only a few arcseconds away 
from each other. On the other hand, the absorption fea- 
tures due to the translucent clouds are very similar in 
the SWAS and HIFI data, in terms of both optical depth 
and width. The angular size of the line of sight clouds 
thus has to be at least comparable to the SWAS beam. 

We estimate the size of the clouds using an indepen- 
dent approach based on the kinematic distance method. 
For a given Galactic longitude and latitude, a given ra- 
dial velocity corresponds to, at most, two heliocentric 
distances. We thus use the Galactic coordinates of the 
corresponding source and the velocity range of each ab- 
sorption component to estimate its size. We convert the 
lower and higher velocities into distances using: 



Vlsr = Ro(0/R - 9 /Ra)sin(l), 



(23) 



where R and Rq are the Galactocentric distances to the 
Sun and the cloud, I is the Galactic longitude of the 
cloud, and 9 and 9o are the Galactic rotation veloci- 
ties at the position of the Sun and the cloud. We use 
a flat rotation curve with 9 = 9q = 254 km s -1 and a 
Gala ctocentric distance to the Sun of 8.4 kpc (|Reid et al.l 
2009). Inherent to the kinematic distance is an ambigu- 
ity between a far and a near distances that both cor- 
respond to the same I and Vlsr- For most PRISMAS 
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lines of sight, there is no ambiguity as the distance to 
the background sources is known. Apart from W49N, 
all the sources are less than 7 kpc from the Sun and the 
translucent clouds are all on the near side. We infer the 
size of each cloud along the line of sight from the differ- 
ence between the kinematic distances that correspond to 
its extreme velocities. For instance, the feature between 
3 and 10 km s _1 towards W51 corresponds to kinematic 
distances between 0.20 and 0.66 kpc. The derived size 
of all the translucent clouds is about a few hundred pc. 
Assuming the projected size on the sky is similar, this 
corresponds to angular scales of a few degrees, much 
larger than the beams of HIFI and SWAS. Naturally, 
this method does not take into account the local veloc- 
ity dispersion of each cloud, therefore the sizes we infer 
should be viewed as upper limits. Additionally, we can- 
not really discuss the dumpiness of the clouds with the 
data in hand. Indeed, a velocity resolution of 200 m/s 
(at 1669 GHz) corresponds to a distance of at least 10 
pc along the PRISMAS lines of sight or to an angular 
size of at least 10 arcminutes, still larger than the HIFI 
and SWAS beams. Nonetheless, if present, those clumps 
or similar short angular scale density variations could 
contribute to the discrepancy between the optical depth 
infered from the 556 and 1669 GHz transitions. 

In this paper, to determine the total column density of 
water, from which we later derive the water abundance 
and the OPR, we therefore decide to use the 1669 GHz 
transition as the best tracer of iVi 01 . 

4.3. On the abundance of water 

The total column density of water for a given gaussian 
component associated with a translucent cloud varies be- 
tween ~ 2 x 10 12 cm -2 and ~ 130 x 10 12 cm' 2 . When we 
integrate over velocity ranges, we usually combine several 
gaussian components. Therefore the total column den- 
sity is usually slightly larger, between ~ 3 x 10 12 cm' 2 
a rid ~ 150 x 1(P cm " 2 . 

iPlume et al.l (|2004T ) analysed the water transitions at 
556 and 547 GHz as observed with the Submillime- 
ter Wave Astronomy Satellite (SWAS) towards W49A, 
which includes W49N. The H 2 spectrum at 556 GHz is 
very similar to the one we present towards W49N except 
that the feature between about and 20 km s _1 , as- 
sociated with the background source, is separated into 
two narrower components in the SWAS observations. 
W e compare in Table El the column densiti es obtained 
bv IPlume et "all (|2TiM )r iGod"ard et al.l ||20H and in the 
present paper, for water and molecular hydrogen, for in- 
terstellar absorption components. We gi ve the values 
we inf er from the gaussian decomposition. IPlume et al.1 
(2004) used CO measurements and assumed a CO abun- 
dance of 10' 4 to infer H 2 column density and water abun- 
dance. We find that there is a rough agree ment, within 
a fact or of a few, between our work and IPlume et al.l 
(2004) regarding the column density of water though the 
SWAS observations seem to indicate the interstellar fea- 
tures are optically thicker than in the HIFI observations. 
The comparison between the H 2 column density derived 
from CO observations and that derived from CH obser- 
vations leads to a larger scattering than for the water 
measurements. There seems to be no systematic trend 
and the ratio between the two estimates of iV(H 2 ) varies 
between 0.16 and 2.3. CH is expected to be a more re- 



liable tracer of the H 2 column than CO in transluscent 
clouds. The use of C O as a tracer of H 2 has been exten- 
sively discussed (e.g. lLiszt et al.ll2010D . 

As a result, th e o-H 2 Q abundance inferred from the 
556 GHz data in IPlume et~a l. (2004) are always larger, 
by a factor of at least a few, than in our work. They find 
abundances between ~ 1 x 10~ 8 and 40 x 10 -8 while we 
obtain abundances between 3.2 x 10 -8 and 7.0 x 10~ 8 . 
Besides the method to derive H 2 abundance, the differ- 
ence in the water vapor column density and abundance 
could arise from the significant difference of the beam 
si ze between SWAS (~ 4') and HIFI (~ 12"). 

INeufeld et all (|2002t ) presented SWAS observations of 
the water transition at 556 GHz towards W51. For 
the feature at about 6 km s _1 , which is not signifi- 
cantly affected by water emission, they infer a water col- 
umn density of 1.9 x 10 13 cm' 2 in good agreement with 

1. 6 x 10 13 cm" 2 we ob tain here. Assuming an OPR of 
3, INeufeld et al.l (|2002t ) claimed the total water column 
density was 2.5 x 10 13 cm" 2 in excellent agreement with 

2.7 x 10 13 cm~ 2 wc infer from the ortho and para ground 
state transitions (from the addition of the 5, 6, and 
7 km s _1 gaussian components in Tabled]). Neufel d et al.l 
(2002) then derived a water abundance of about 10~ 8 
relative to Hi while we derive an abundance of 6 x 10 ~ 8 
relative to molecular hydrogen. Those water abundances 
are in agre ement with each othe r within less than a fac- 
tor 2 since IGodard et a l. (2012) report for this feature 
column densities of atomic hydrogen (14.8 x 10 20 
and molecular hydrogen (3.5 — 4.3 x 10 20 cm 
respond to a Hi-to-H 2 ratio of 3.4 — 4.2. 

We use the H 2 column densities that IGodard et al.l 
(2012) obtained from CH absorption line, assuming 
HF/CH = 0.4 and HF/H 2 =3.6 x 10~ 8 , since the HF ab- 
sorption line are more often saturated. We r eport these 
H 2 column densities values in Table Figure 9(a) shows 
the total H 2 column density as a function of the H 2 col- 
umn density for nine translucent clouds along the PRIS- 
MAS lines of sight. We select those as (1) the velocity 
position and width of the water emission fr om the source 
is not an obvious contamination and (2) IGodard et al.l 
(|2012f) provide a column density for H 2 (all the sources 
except W28A). The nine clouds span a range of column 



Table 5 

Comparison with the results from Plume ct al. (2004) towards 
W49 

Absorption features 
(km s — 1 ) 



cm ) 
~ 2 ) that cor- 



33.5 39.5 


53.5 


59.6 


63.3 68.0 


jV(o-H 2 0) 


(a) 


32 >55 


30 


>100 


>25 13 


AT(o-H 2 0) 


00 


53 


26 


39 


34 


Ratio 




>1.6 


1.2 


>2.6 


>1.1 


AT(H 2 ) 


(a) 


0.8 15 


1.3 


4.6 


28 1.6 


AT(H 2 ) 


(c) 


5.9 6.6 


3.7 




22.6 


Ratio 




0.16 2.3 


0.35 




1.5 


X(o-H 2 0) 


(a) 


40 >3.6 


23 


>22 


>0.92 8.1 


X(o-H 2 0) 


(b) 


4.2 


7.0 




3.2 



Note. — To remain consistent with IPlume ct al. (20041), we 
assume here that N(o— H 2 0) = Ni 01 and that Ni 01 is given by 
the transitions at 556 GHz. The water column density is given 
in units of 10 12 cm -2 , that of H 2 in units of 10 20 cm -2 , and 
the relative abunda nce of water in un its of 10~ 8 . R e ferenc es: (a) 
IPlume et all 12001 1. (b) this work, (c) IGodard et"all (1201211 



Water absorption in the Milky Way 



17 




10 15 20 
N(H 2 ) /10 20 (cm" 2 ) 

(a) The red solid line is the best linear fit to the data. 



30 




40 60 80 100 120 
N H /10 20 (cm" 2 ) 

(b) The red solid line is the best linear fit to the data, 
show the features for which we only have a lower limit 



Arrows 
for N H 



Figure 9. Total column density of water as a function of the column density of molecular hydrogen and hydrogen nuclei. 
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Figure 10. Abundance of water relative to molecular hydrogen and hydrogen nuclei as a function of 7V(H2) or Njj. 
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Figure 11. Abundance of water relative to molecular hydrogen and hydrogen nuclei as a function of the molecular fraction of hydrogen. 

density that is similar to our whole sample, from about 10 x 10 12 cm -2 to 120 x 10 12 cm -2 for H2O. There is 



2 x 10 2U cm" 2 to 30 x 10 2U cm -2 for H 2 , and from about 



a clear linear correlation between the column densities 
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which indicates a remarkably constant abundance of H2O 
from low to high H2 column density clouds. A linear fit 
to the data gives a slope of (5.0±0.4) x 10~ 8 if we assume 
there are no water molecules in the excited levels and the 
uncertainty is thus only given by the uncertainty on the 
ground state column densities. We also show the total 
column density of water as a fu nction of the t o tal h ydro- 
gen column density Nh from iGodard et al.l (|2012t ) for 
the same nine clouds, though we only h ave lower-limits 
of Nh for two of them (see Figure 9(b) ). Again, we find 
a clear correlation, though not as good as with N(H 2 ) 
suggesting that water is a better tracer of H2 than hy- 
drogen. A linear fit gives a slope of (1.4 ± 0.2) x 10~ 8 if 
we use the uncer t ainty on the water ground levels only. 
iHollenbach et all (|2012f) modeled the production of oxy- 
gen bearing molecules (OH + , H20 + , and H30 + ) in inter- 
stellar clouds and found that, for a density of 100 cm -3 , 
at a depth between Ay ~ 0.1 and Ay ~ 1 mag, the 
abundance of water relative to hydrogen nuclei was al- 
most constant at ~ lx 10 -8 . To obtain this abundance, 
IHollenbach et all ()2012h used an average incident diffuse 
interstellar radiation field (x = 1) and a high cosmic ray 
ionization rate of 2 x 10 -16 s" 1 . Our measurements of 
the water abundance relative to hydrogen nuclei are in 
very good agreement with their results. 

We then explore possible variations of the water abun- 
dance, whose values are indic ated in Ta ble [6] for each 
cloud. First, we show in Figures 10(a) and 10(b) the vari- 
ations of the water abundance, relative to molecular hy- 
dogen, or hydrogen nuclei, as a function of N(H 2 ), or Nh- 
Neither plot show any sign of a correlation. The average 
value for JV(H20)/AT(H2), taking into account the only 
uncertainty of each data point due to the ground state 
column densities, is (4.8 ±0.3) x 10" 8 . The average value 
for N(R 2 0)/N H is (0.8±0.1) x 10~ 8 whatev er uncertain- 
ties we use. Then, we show in Figures 11(a) and 11(b) the 
abundance of water relative to molecular hydrogen, or 
hydrogen nuclei, as a function of the molecular fraction of 
hydrogen (/(H 2 ) = 2N(K 2 )/N H ). If N(K 2 0) was dire- 
city proportional to Nh , then we should see a correlation 
between N(R 2 0)/N{R 2 ) and /(H 2 ) = {N H /2N{R 2 ))- 1 . 
Instead, we see that N (H2O) / 'N '(Ha) does not exhiti 
any trend as a function of /(H2). On the other hand, 
N(R 2 0)/N H increases as /(H 2 ) = 2N(R 2 )/N H in- 
creases, which confirms that water is a better tracer of 
H2 than it is of hydrogen. The water abundance relative 
to hydrogen nuclei increases from less than 0.5 x 10 ~ 8 
to about 1.5 x 10 ~ 8 as /(H2), which is th e range of val- 
ues predicted by IHollenbach et all ()2012h as the visual 
extinction, therefore the molecular fraction of hydrogen, 
increases from less than 0.1 mag to about 1 mag. This 
confirms that N(R 2 0) is a good tracer of N(H 2 ) in the 
translucent section of clouds. 

We finally explore the possiblity of a variation of the 
water abundance with Galactocentric distance as the 
oxygen abundance in the Galaxy has been shown to 
decrease a few 0.01 dex/kpc with Galact ocentric dis- 
tance over significantly larger range (e.g. iBalser et al.l 
[20lH iHenrv et al.ll201ll ). Therefore, we associate each 
cloud with its Galactocentric distance, inferred from the 
kinematic distance method, to determine whether the 
abundance of water also follows a gradient through the 
Galaxy. We use the whole range of velocity covered by 



Table 6 

Abundance of water relative to molecular hydrogen and to 
hydrogen nuclei. 



Cloud 



jV(H2 0)/Ar(H 2 ) N(H 2 Q)/N H 



DR21(OH) (4 to 18 km s" 1 ) 


' 0.4 


<0.77 


G34.3+0.1 (8 to 16 km s" 1 ) 


7.7+I 7 


1 0+O.3 
^•^0.1 


G34.3+0.1 (22 to 34 km s" 1 ) 


4 2+ 2 ' 2 


041 +0.17 


G34.3+0.1 (42 to 55 km s" 1 ) 


<10.4 




W33(A) (21 to 26 km s" 1 ) 


4+ 2 


n oo+0. 16 


W33(A) (26 to 31.5 km s" 1 ) 


4.8 +2 8 3 


i-2 + ;- 4 


W49(N) (26 to 48 km s" 1 ) 


5+ 2 


o.8 + r 


W49(N) (48 to 73 km s" 1 ) 


4 -!+ - 9 


1 0+ 01 


W51 (3 to 15 km s" 1 ) 


6+ 3 


94+ 1 


W51 (44 to 51 km s" 1 ) 


12+ 15 


<1.1 



Note. — JV(H 2 0) are in units of Id 1 



JV(H 2 ) and N H 



are in units of 10 cm 



Abundances are in units of 10 



each absorption feature to determine its Galactocentric 
distance and estimate the uncertainty on this distance. 
We find that the nine translucent clouds for which we 
measure a wate r abun dance are located between 5.5 and 
8.5 kpc. Figure 12(a) shows that neither the H2O abun- 
dance relative to molecular hydroge or that relative to 
hydrogen nuclei varies significantly with the distance to 
the Galactic center. Only additional measurements over 
a wider range of Galactocentric distances might help bet- 
ter constrain a trend, if any. 

4.4. The OPR as a tracer of the gas history 

Figure [TBI shows the OPR for each interstellar gaussian 
absorption feature in the PRISMAS spectra as a function 
of the H2O column density. No obvious trend can be de- 
tected and a linear fit to the data leads to a slope of less 
than 10 -2 . The weighted average of the OPR is about 
2.8±0.1. Most features are in close or very close agree- 
ment with the high temperature limiting value of 3. In- 
deed, for 34 out of 41 interstellar absorption components, 
the observed OPR is within 3-er of a value of 3.0. Only 
a few gaussian components are either significantly below 
or above 3. Among them, the most extreme values corre- 
spond to the features at 55 km s _1 towards W49N, and 
at 29 km s _1 towards W33A. These two features have in 
common that they are blended with other stronger ab- 
sorption features, including features associated with the 
background source. As a result, the column density and 
OPR of water derived for those absorption components 
might only be correct when summed over a larger veloc- 
ity range that includes all the blended features. 

We then analyse the OPR derived from the column 
density curves. Hereafter, we use OPR 9 as the estima- 
tor of the OPR, rather than OPR + , since we expect the 
water column density in the excited levels to be negligi- 
ble, assuming T ex ~ 5 K. In total, we measure the water 
OPR for 13 translucent clouds. For these 13 clouds, the 
average OPR is 2.9±0.1, in very close agreement with 
the high temperature limiting value of 3. Out of the 
13 clouds, 10 have an OPR less than 3-cr away from a 
value of 3. One of the three other clouds has an OPR 
value above 3 and correspond to the narrow absorption 
component at ~ 44 km s _1 towards W33(A). The OPR 
inferred for the remaining two clouds is 2.3±0.1, 2.4±0.2. 
These are the clouds towards W49N at ~ 40 km s _1 and 
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Figure 12. Abundance of water relative to molecular hydrogen and hydrogen nuclei as a function of the Galactocentric distance. Each 
data point represents one absorption feature. The label indicate the background source and the velocity for each absorption feature. 

ponents to that obtained for each cloud, combining the 
gaussian components accordingly. We find there is al- 
most a perfect agreement between the result from the two 
methods (see Figure[T4|). Some discrepancy is notable for 
a couple of clouds where the gaussians give a lower than 
3 OPR while the column density curves leads to an OPR 
in agreement with a value of 3. The two clouds towards 
W49N are the only ones that have an OPR significantly 
lower than 3 for both methods. This statistical analy- 
sis on 13 cloud s confirms and extends the results from 
iLis et al.1 (|2010f ) who analyzed two PRISMAS sightlines 
and found some values of the water OPR in agreement 
with the high temperature limiting value of 3 and some 
values that were significantly below this limit. 

Figure [15] shows the OPR derived for each translucent 
cloud using the column density curves as a function of 
the total water column density, the molecular hydrogen 
column density, the hydrogen nuclei column density, and 
the hydrogen molecular fraction. We assume the wa- 
ter column density inferred for the excited levels is only 
due to the noise. There seems to be a slight decrease of 
the OPR as the column density increases, whether it is 
that of H2O H2 or hydrogen nuclei. If true, this trend 
should imply a decreas e of the OPR as /(H2) increases. 
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Figure 13. Water ortho-to-para ratio as a function of the total 
water column density both derived from the gaussian absorption 
components. The dashed blue line shows the high temperature 
limiting value of 3. 
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Figure 14. Comparison between the water OPR derived from 
the gaussian components and that derived from the column density 
curves, assuming there are no water molecules in the excited levels. 
The dashed blue line shows where the two OPR ratios are equal. 

- 60 km s _1 . 

We compare the OPR derived from the gaussian corn- 



ways appear on the high end of the column density axis. 
However, the limited number of clouds and amplitude of 
variations though prevent us from drawing a definitive 
conclusion regarding a correlation between a low OPR 
and a high column density. 

The two translucent clouds for which we measure an 
OPR of ^2.3 are observed towards W49N, at the veloc- 
ities of ~ 40 and ~ 60 km s _1 (see Figure |2"0|) . These 
clouds span a large velocity range of about 20 km s" 1 
and combine many velocity sub-components, which ex- 
plains why we derive a large total column density of wa- 
ter. Neither of them show evidence of saturation though 
the peak of the sub-component at 40 km s^ 1 seems to 
just reach the zero level of intensity at 556, 1113, and 
1669 GHz. However, none of the components exhibit a 
plateau at the zero level of intensity. Additionally, the 
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Figure 15. Water ortho-to-para ratio as a function of the water column density, the molecular hydrogen column density, the hydrogen 
nuclei column density, and the hydrogen molecular fraction. The dashed blue line shows the high temperature limiting value of 3. 



background emission seems well constrained over the ve- 
locity range of these two clouds as the emission features 
from the background source peak at ^5 km s~ x and ap- 
pear to be limited to velocity smaller than ~ 30 km s" 1 . 
One peculiarity of these clouds is due to the position of 
the background source, W49N, which is the only PRIS- 
MAS source for which there might be an ambiguity be- 
tween the near and the far distances derived from the 
kinematic distance method. Figure Q] shows that the ve- 
locities of 40 and 60 km s _1 can either be related to 
structures at about 3 kpc or 9 kpc from the Sun. The 
Perseus and Sagittarius arms are located at those ap- 
proximate distances towards W49N. Consequently, the 
absorption features we see towards W49N at those ve- 
locities might be a combination of near and far distance 
clouds. Nonetheless, the inferred OPR values is signif- 
icantly below 3 at least in some parts of these clouds. 
The cloud at ~ 40 km s _1 is clearly detected in the ob- 
servations of nitrogen hydrides, es pecially in NIT, fo r 
which the OPR is peculiar as well (jPersson etaDESil- 
At millimeter wavelengths, the cloud is also detected in 
H 13 CO + , which shows that it has a high column density 
and properties in the high end of the transluscent cloud 
regime approaching dark clouds/dense cores. The situa- 



tion is less clear for the feature at ~ 60 km s _1 which is 
probably a mixture of several environments. 

We relate the low OPR values to equilibrium tempera- 
tures below 50 K. An OPR of ^2.3 corresponds to water 
vapor in thermal equilibrium at ^25 K. The low equi- 
librium temperature could be a result of the formation 
of the water molecules in cold gas, on cold dust grains, 
or the thermalisation with cold gas of water molecules 
formed at higher temperature. When water forms in 
the cold gas phase through H 3 + recombination, the 
large amount of energy released in this exothermic re- 
action should lead to spin equilibration of the molecules. 
Similarly, when water molecules are desorbed from dust 
grain ice mantles, the energy required to do so is suffi- 
cient to populate many ort ho and para levels and lead to 
an OPR of 3. Nonetheless. iLimbach et ail (|2006f) showed 
that the OPR of water vapor released from a cold surface 
(e.g. ice mantle of dust grain) could reflect the temper- 
ature of the surface and of the gas phase molecules close 
to the surface. Indeed, the energy required to desorb 
water molecules from dust particles can be dissipated 
into th e grain as well as the molecule. iHollenbach et all 
(2009) also showed than in diffuse regions, molecules of 
water could be released in a frozen state with a spin 
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temperature equal to that of the grain, which is smaller 
than that of the gas in diffuse and translucent clouds. 
Therefore, the released molecules might carry the in- 
formation about the temperature where they have been 
formed. However, photodesorption may be the only ef- 
ficient way to remove H2O molecules from cold dust 
grains and, depending of the exact mechanism (either 
H2O dissociation followed by recombination and escape, 
or H2O dissociation followed by ejection of neighbor- 
ing H2 O molecule), the OPR may or m ay not be pre- 
served (jAndersson fc van Dis hoeck 2008). Additionally, 
at high temperature (e.g. hot cores, outflows), the gas 
phase reactions 0(H2, H)OH(H2, H)H20 can also from 
water molecules that are expected to have an OPR of 
3. Whether it was formed in cold gas, warm gas or des- 
orbed from a dust grain ice mantle, a water molecule in 
the gas phase will subsequently thermalize through colli- 
sions with atoms and molecules. During these collisions, 
proton exchanges may occur and the OPR of the water 
molecules will then carry information about the temper- 
ature of the medium with which they last interacted. 
Whether the observed OPR of water molecules relate to 
their formation or their more recent history depends on 
the efficiency of the thermalisation and spin conversion 
processes. 

In order to allow such spin conversion, a sufficient 
number of collisions with ions are required. This may 
relate to a large density or a long timescale. For in- 
stance, within a 10 cm -3 gas, protonated ions with 
a fractional abundance of 10 -8 , and a rate coefficient 
of 10" 9 cm _3 s _1 , the OPR equilibriation time scale is 
about 3 x 10 5 years. In a translucent cloud with a 
characteristic density of 100 cm" 3 , and other param- 
eters as above, the OPR equilibriation time scale is 
about 3 x 10 7 years. The low OPR values we detect 
in the two clouds towards W49N implies that the tem- 
perature of the cold gas (~ 25 K) is below the ex- 
pect ed temperature in tr anslucent clouds (~ 50 K, see 
e.g. iRachford et al.l 120021) . We characterize the clouds 
along the lines of sight as translucent because of the 
extinction (1-5 ma g) derived from the hydrogen col- 
umn densities that iGodard et al.1 (|2012l ) have reported. 
In the two clouds towards W49N, the total hydro- 
gen column densities are (86.7±10.6) x 10 12 cm -2 and 
(116.6±10.2)x 10 12 cm" 2 . The corresponding extinctions 
are of 4.6±0.6 and 6.2±0.5 mag using R y = 3.1 and 
N H /E(B-V) = 5.8 x 10 21 H cm^mag" 1 (|Bohlin et al.l 
[197a IRachford et al.ll200l . These are sufficient for wa- 
ter molecules to be stuck on dust grains as the extinc- 
tion threshold to for m water ice mantles has been shown 
to be about 3 mag (fWhittet et all 120011 ). Therefore the 
water molecules released from the grains can thermalize 
with the cold gas in a few 10 5 years and have an OPR 
lower than 3. If those water molecules later move to the 
translucent regions of the clouds, they will not thermalize 
in less than a few 10 7 years. This indicates that the wa- 
ter molecules for which we measure an OPR significantly 
below 3 have (1) either thermalizcd with the cold gas 
or been released cold from grains within denser, colder 
regions of the clouds and (2) have not yet thermalized 
at the temperature of the warmer, translucent regions of 
the clouds. 

In the present study, the low value of the OPR found 



in several translucent clouds may be related to any of the 
aforementioned interpretations, either based on the for- 
mation mechanisms or the thermalization with gas colder 
than ~50 K. The conclusion is nevertheless that the wa- 
ter molecules have been in a cold environment, be it in 
the gas phase or on the dust grains. 

5. CONCLUSION 

For the first time, Herschel/HIFI allows the observa- 
tion of the three ground state transitions of H2O and 
H2 8 as well as the first three excited transitions of 
H2O in absorption towards six continuum sources in the 
Galactic plane. The continuum sources are detected at 
all frequencies in self- absorption associated to an enve- 
lope or outflow. They are also detected in broad water 
line emission at almost all frequencies. We detect the 
translucent clouds associated with all the Galactic arms 
along those lines of sight only in the ground state tran- 
sitions of both ortho and para-H^O. 

For each transition and sightline, we model the con- 
tinuum emission, taking into account the contribution of 
the water lines and decompose the optical depth into a 
sum of gaussian profiles. From there, we infer the col- 
umn density of water in the ground state levels of H2O, 
and upper limits for the excited levels of H2O and the 
ground state levels of H^O, for each interstellar absorp- 
tion feature. The total column density of water ranges 
from a few 10 12 to about 10 14 cm -2 . We compare the 
column density of water to that of molecular hydrogen 
and find that the water abundance is remarkably con- 
stant at ~ 5 x 10" 8 . Relative to hydrogen nuclei, the 
water abundance is ~ 10~ 8 , a value predicted by mod- 
els with high cosmic ray ionization rate. No correlation 
with the hydrogen molecular fraction or Galactocentric 
distance are found. Water therefore seems to be a good 
tracer of H2 in translucent clouds. 

For the first time, the excited transitions of H2O are 
observed. Though no translucent clouds are detected in 
absorption in those transitions, we infer upper limits on 
the column density of H2O in the corresponding lower 
levels. We then derive the lower limits on the fractional 
abundance of water molecules in the ground states and 
find it is usually above 85%. Using this parameter, we set 
constraints on the excitation temperature (T ex < 12 K). 

For the first time too, both ground state transitions of 
ortho water have been observed. Some discrepancies are 
revealed and we suggest several interpretations. First, 
assuming that the excitation temperature is low (T ex -C 
hv/k) may be wrong. We find that with T ex ~ 5 K, 
both ortho ground state transitions are consistent, and 
that the water fractional abundance in the excited levels 
is significantly lower than 1%. These values are both 
in agreement with the upper limits derived from the 
observations of the excited transitions. The excitation 
temperature then enables us to provide an upper limit 
for the density in the clouds (< 10 4 cm" 3 , which is in 
agreement with the typical density in translucent clouds 
(~ 100 cm -3 ) but does not rule out the existence of more 
extinguished regions in the clouds. Other interpretation 
of the discrepancy between the two ortho ground state 
transitions may arise from the significant difference in 
the continuum strength at both frequencies, and/or the 
presence of clumps in the translucent clouds whose char- 
acteristic size would imply different filling factor within 
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beams of 14.5 and 41" though a comparison with SWAS 
observations makes this last possibility unlikely. 

We finally derive the water OPR for each interstel- 
lar absorption component. Most clouds have a ratio in 
agreement with the high temperature limiting value of 3. 
However, we also found that a few features correspond to 
an OPR significantly below 3. We discuss several inter- 
pretations, in view of the properties of these clouds which 
suggest the presence of high extinction regions (~ 5 mag) 
where the water could be stuck on cold grains. Whether 
the water molecules were released with a "cold" OPR or 
quickly thcrmalizcd with the cold gas in the dark sections 
of the clouds remains unclear though we claim they have 
not yet thermalized at the temperature of the translucent 
section of the clouds. 

An analysis of more subtle variations within each ab- 
sorption features, allowed by the sensitivity and resolu- 
tion of HIFI on Herschel, correlated with the many other 
molecular tracers of the ISM properties that have been 
observed within the PRISMAS Guaranteed Time Key 
Program may eventually allow us to discriminate among 
the proposed explanations of a low OPR in relatively 
warm regions. 
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APPENDIX 
OBSERVATIONS 



In this section of the appendix, we gather all the spectra (nine transitions) for all the PRISMAS sources except 
W51, which is presented in the main part of the paper. 
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Figure 16. Same as Figure[3]but towards DR21(OH). 
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Figure 17. Same as Figure|3]but towards G34. 3+0.1. 
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Figure 18. Same as Figure|3]but towards W28(A). 
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Figure 19. Same as Figure |3]but towards W33(A). 
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Figure 20. Same as Figure|3]but towards W49(N). 
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RESULTS 

In this section of the appendix, we give the column densities of H2O, H^O, and H2 for each absorption feature, the 
inferred abundance, and OPR., for all the PRISMAS sources except W51, which is presented in the main part of the 
paper. 



Table 7 

Same as Table|3]but towards DR21(OH). 
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Same as Tablc[3]but towards G34.3+0.1. 
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60.6 


5.62 


3.49 


47 ib 2 


IfifiQ 


11.7 


1.98 


2.53 


24 zb 1 




1 A 1 


i in 
± . iy 


n ^oa 
u.ou^t 


9 ft 4- f) 9 

z.o m u.z 




27.4 


1.15 


1.66 


9.1 ± 0.1 




27.9 


5.46 


0.251 


6.7 ± 0.3 




45.5 


4.56 


0.291 


8 ± 1 




48.7 


2.16 


1.22 


11 ± 3 




52.6 


3.20 


2.46 


38 ± 2 




60.7 


6.58 


3.57 


110 ± 4 


547 GHz 


60.6 


2.68 


0.124 


3 ± 1 


1101 GHz 


61.4 


2.97 


0.201 


1.4 ± 0.2 


1655 GHz 


61.6 


10.9 


0.255 


12 ± 2 



Table 9 

Same as Table [3] but towards W28(A). 



Transition 


"0 


FWHM 


To 


AT;(H 2 0) 


(GHz) 


(km s" 1 ) 


(km s _1 ) 




(xlO 12 cm" 2 ) 


987 GHz 


-1.25 


12.4 


0.377 


43 ± 2 


1661 GHz 


-1.05 


12.8 


0.695 


123 ± 2 


556 GHz 


-1.35 


19.8 


0.394 


35.6 ± 0.4 




1.29 


9.40 


0.575 


25 ± 2 




5.97 


3.82 


0.990 


17.5 ± 0.6 




12.1 


2.54 


2.77 


32.7 ± 0.1 




14.4 


1.23 


1.01 


5.8 ± 0.3 




19.7 


3.08 


5.92 


85 ± 1 




23.8 


2.83 


0.840 


11.1 ± 0.2 


1113 GHz 


-1.36 


36.3 


0.469 


40.9 ± 0.1 




1.20 


10.0 


0.809 


19 ± 1 




6.12 


3.49 


0.896 


7.3 ± 0.6 




12.1 


2.46 


1.87 


10.7 ± 0.2 




14.5 


0.983 


0.915 


2.1 ± 0.1 




19.6 


2.86 


4.42 


30 ± 1 




23.6 


3.41 


0.627 


4.9 ± 0.4 


1669 GHz 


-1.43 


38.1 


0.556 


94 ± 5 




1.22 


9.84 


0.971 


48 ± 5 




6.05 


3.92 


1.37 


25.2 ± 0.5 




12.2 


2.52 


3.06 


37 ± 1 




14.7 


0.885 


1.39 


5.8 ± 0.1 




19.7 


2.77 


7.94 


103 ± 3 




23.7 


3.17 


0.897 


13.4 ± 0.3 


547 GHz 


-1.27 


23.9 


0.147 


17 ± 1 


1101 GHz 


-2.11 


23.9 


0.149 


8.2 ± 0.9 


1655 GHz 


-1.79 


23.9 


0.379 


42.0 ± 0.1 
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Table 10 

Same as Table |3] but towards W33(A). 



Transition 


■"0 


r Wrilvl 


TO 


AT /TU f\\ 

i\l(i±2U) 


[\jjtiZ ) 


S J 


S J 




^XiU cm J 




O4.0 


o.yu 


U.z / O 


1U.U ± U.o 


-I /;t?1 /"I T_J 

lool GHz 


33.6 


6.96 


0.392 


46 ± 7 


556 GHz 


24.0 


1.61 


0.698 


5.3 ± 0.3 




29.5 


3.92 


1.28 


22 ± 5 




32.5 


4.31 


0.882 


17.2 ± 0.5 




33.7 


2.17 


3.97 


40 ± 1 




37.0 


2.14 


4.10 


40.9 ± 0.8 




39.2 


1.21 


6.26 


35 ± 3 




40.9 


1.97 


1.03 


8 ± 3 




44.2 


0.821 


0.218 


0.8 ± 0.1 


1113 GHz 


24.1 


1.69 


0.512 


2.0 ± 0.1 




29.6 


3.93 


1.24 


11.5 ± 0.9 




33.5 


1.34 


2.08 


6 ± 1 




34.4 


4.78 


2.10 


23 ± 2 




37.1 


1.20 


3.68 


10.3 ± 0.9 




39.0 


1.74 


4.04 


16.3 ± 0.3 




41.4 


1.41 


0.703 


2.4 i 0.2 




A A 1 

44. z 


0.731 


0.295 


0.5 i 0.1 


IDOy LrrlZ 


Z4.o 


z.Ul 


U.O /O 


4 ± 1 




29.3 


2.48 


1.73 


ly.o m u.o 




31.0 


6.45 


0.960 


21.8 ± 0.3 




33.9 


3.54 


3.26 


58 ± 2 




37.5 


2.14 


3.39 


33 ± 6 




39.4 


1.87 


2.88 


20 ± 13 




41.6 


1.75 


0.996 


8 ± 2 




44.3 


0.827 


0.537 


2.2 ± 0.1 


547 GHz 


33.1 


4.40 


0.135 


3.1 ± 0.1 


1101 GHz 


34.3 


2.62 


0.0552 


0.7 ± 0.7 


1655 GHz 


33.3 


7.00 


0.208 


7.4 ± 0.2 



Water absorption in the Milky Way 



Table 11 

Same as Table |3] but towards W49(N). 



Transition 


Vq 


r WrllVl 


TO 


AT /u r\\ 


\\jrVL ) 


(km s — ) 


(km s — ^) 




(Xiu cm ) 


9o7 GHz 


11.9 


7.78 


0.398 


26 ± 4 


1661 GHz 


9.39 


10.4 


1.22 


175 ± 2 


556 GHz 


3.71 


1.87 


0.195 


1.9 ± 0.4 




6.59 


2.15 


1.34 


13.5 ± 0.4 




11.5 


9.13 


1.84 


71 ± 9 




16.1 


3.08 


5.06 


88 ± 5 




33.6 


2.54 


1.02 


11.8 ± 0.5 




39.2 


5.93 


0.692 


20.3 ± 0.5 




39.4 


1.49 


3.08 


21 ± 1 




51.3 


1.76 


0.279 


2.3 ± 0.1 




55.1 


4.71 


1.07 


23.4 ± 0.7 




59.8 


2.88 


2.86 


39 ± 2 




62.8 


2.66 


2.21 


27.4 ± 0.9 




65.8 


1.66 


0.221 


1.7 ± 0.2 




68.3 


3.62 


0.275 


4.6 ± 0.3 


1113 GHz 


6.22 


1.07 


0.538 


2.0 ± 0.3 




10.6 


10.8 


2.66 


68 ± 3 




15.6 


2.75 


9.51 


61 ± 30 




33.5 


1.67 


0.675 


2.6 ± 0.2 




38.3 


9.81 


0.405 


9.3 ± 0.7 




39.4 


1.47 


2.62 


9.1 ± 0.8 




51.2 


1.29 


0.194 


0.6 ± 0.1 




55.4 


6.30 


0.675 


9.9 ± 0.3 




60.1 


3.19 


1.73 


12 ± 1 




62.9 


1.89 


1.36 


5.9 ± 0.5 




65.2 


2.42 


0.182 


1.0 ± 0.7 




68.3 


4.52 


0.202 


2.2 ± 0.2 


1669 GHz 


6.70 


1.39 


0.764 


6 ± 3 




10.5 


11.8 


2.53 


138 ± 5 




16.1 


2.53 


17.5 


202 ± 47 




33.8 


1.55 


0.916 


6.7 ± 0.1 




Ol .0 


1 1 o 

11. Z 


U.40y 


ZO.O ± U.Z 




39.5 


1.84 


2.15 


i c e. _|_ n q 




51.2 


1.42 


0.308 


2.1 ± 0.1 




55.0 


4.75 


0.836 


18.5 ± 0.1 




60.3 


4.00 


1.99 


37 ± 1 




63.2 


1.90 


1.47 


13.1 ± 0.3 




65.7 


2.02 


0.255 


2.5 ± 0.3 




68.6 


3.49 


0.273 


4.4 ± 0.1 


547 GHz 


12.0 


7.26 


0.222 


7.1 ± 0.5 


1101 GHz 


10.5 


9.68 


0.126 


3.1 ± 0.6 


1655 GHz 


10.3 


12.0 


0.357 


20.0 ± 0.5 



Table 12 

Same as Table|4]but towards DR21(OH). 



V 

(km s" 1 ) 


AT(H 2 0) 


OPR 


V 


JV(H 2 0) 


/s(H 2 0) 


OPR+ 


OPR 3 


AT(H 2 ) 


X(H 2 0) AT(H^O) 


8 


34 ± 2 


4.1 ± 0.8 
















9 

12 


33.5 ± 0.5 
5.7 ± 0.1 


3.0 ± 0.1 
2.9 ± 0.1 


4 to 18 


72+1 


>0.93 


2.8±0.5 


2.9±0.3 


15.2±1.0 


4.7i°-j <2.5 


15 


5.1 ± 0.1 


2.8 ± 0.1 

















Note. — The water column densities 7V(H 2 0) are in 10 12 cm 2 . 
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Table 13 

Same as Tablegjbut towards G34.3+0.1. 



V 

(km s" 1 ) 


^(HaO) 


OPR 


V 


iV(H 2 0) 


/ S (H 2 0) 


OPR+ 


OPR 9 


AT(H 2 ) 


X(H 2 0) 


N(Hl»0) 


12 
14 


32 ± 1 
3.7 ± 0.2 


3.1 ± 0.2 

3.2 ± 0.2 


8 to 16 


36 1 +s 


>0.84 


2.8±0.8 


3.0±0.2 


4.7±0.6 


7 7+2.7 
'•'-1.1 


<0.3 


27 
28 


12.3 ± 0.1 
9.2 ± 0.4 


2.8 ± 0.1 
2.6 ± 0.4 


22 to 34 


22 5+ 9 


>0.73 


3±1 


2.9±0.3 


5.3±0.6 


4 2+ 2 - 2 


<0.7 


46 
49 
52 


11.2 ± 1.3 

14 ± 3 
51 ± 2 


2.4 ± 0.4 

2.8 ± 0.8 

2.9 ± 0.2 


42 to 55 


78^2° >0-73 
Group of absorption 


3.2±0.9 2.9±0.1 >10.4 <10.4 
features partly contaminated by the source. 


<3.4 



Note. — The water column densities jV(H2 0) are in 10 12 cm 2 . 



Table 14 

Same as Table [4] but towards W28(A). 



V 

(km s~ 




Af(H 2 0) 


OPR 


V 


7V(H 2 0) / 9 (H z O) 


OPR+ OPRf 7V(H 2 ) 


X(H 2 0) 


N(HI«0) 


1.0 




67 ± 5 


2.6 ± 0.4 




Group of absorption 


features contaminated by the 


source. 




6.1 




32.6 ± 0.8 


3.4 ± 0.3 












12 
15 


47 ± 1 
7.9 ± 0.1 


3.4 ± 0.2 
2.8 ± 0.1 


9.5 to 16 


73+ 7 2 >0.93 


3.0±0.4 3.0±0.2 




<4.7 


20 




133 ± 3 


3.4 ± 0.2 




14ll^3 >0.98 








24 




18.3 ± 0.5 


2.7 ± 0.3 


16 to 28 


3.3±0.6 3.4±0.5 




<2.4 



Note. — The water column densities 7V(H 2 0) are in 10 12 cm 



Table 15 

Same as Table g] but towards W33(A). 



V 

(km s" 1 ) 


7V(H 2 0) 


OPR 


V 


7V(H 2 0) / 9 (H 2 0) 


OPR+ OPRs iV(H 2 ) 


X(H 2 0) 


iVtHfO) 


24 


6.3 ± 1.3 


2.1 ± 0.7 


21 to 26 


8.7^02 >0.74 


3±1 2.9±0.1 2.4±0.8 


3.eHi% 


<1.2 


29 


31 ± 1 


1.7 ± 0.2 


26 to 31.5 


44.3+i° 6 >0.75 


3.3±0.9 3.0±0.2 9.2±1.2 


4 8 +2 - 3 

4 -°-0.8 


<2.5 


37 


44 ± 6 


3.2 ± 0.9 












39 


37 ± 14 


1.3 ± 0.9 




Group of absorption 


features contaminated by the 


source. 




41 


10 ± 2 


3.4 ± 1.0 












44 


2.7 ± 0.1 


4.3 ± 0.2 


43.2 to 45.5 


2.9+o.t >0.90 


4.2±0.8 4.3±0.3 




<0.2 



Note. — The water column densities 7V(H 2 0) are in 10 12 cm 2 . 



Table 16 

Same as Table H] but towards W49(N). 



V 

(km s" 1 ) 


N(H 2 0) 


OPR 


V 


JV(H 2 0) / 9 (H 2 0) 


OPR+ OPR» AT(H 2 ) 


X(H 2 0) 


AT(H^O) 


7 


8 ± 3 


3 ± 2 












16 


263 ± 56 


3 ± 2 




Group of absorpt 


on features contaminated by the 


source. 




34 


9.3 ± 0.2 


2.6 ± 0.2 




73.3lg 9 > 0.94 








38 


33.0 ± 0.7 


2.6 ± 0.2 


26 to 48 


2.2±0.3 2.3±0.1 13.6±3.1 


<1.9 


39 


27.6 ± 0.8 


2.0 ± 0.2 












51 


2.7 ± 0.1 


3.2 ± 0.2 












55 


28.3 ± 0.3 


1.9 ± 0.1 












60 


50 ± 1 


2.9 ± 0.3 


48 to 73 


112l| > 0.96 


2.3±0.3 2.4±0.2 26.3±3.8 


4.3ig : » 


<1.8 


63 


19.0 ± 0.6 


2.2 ± 0.2 










65 


3.5 ± 0.8 


3 ± 2 












68 


6.6 ± 0.2 


2.0 ± 2.2 













Note. — The water column densities Af(H 2 0) are in 10 12 cm 2 . 



